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ABSTRACT 



A three-dimensional analysis of the intensity distribution of backward optical 
transition radiation has been performed. The effects of variations in electron energy 
and beam divergence and on material properties such as dielectric permittivities and 
the resultant coherence length upon the angular distribution and polarization of 
optical transition radiation has been investigated. A surprising obeservation 
important to the use of optical transition radiation as a diagnostic tool for high 
energy electron beams is the behavior of the perpendicular component of the 
intensity. In contrast to low energies where the parallel component dominates, at 
electron energies above 200 MeV, the perpendicular component dominates. This 
requires the use of a polarization filter to diagnose particle beam properties at high 
energies. 
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I. INTRODUCTION 



Transition radiation is produced whenever a uniformly moving charged 
particle passes from one medium into another. In introducing the concept in 1944, 
Frank and Ginsburg [Ref. 1] noted that the intensity, polarization and angular 
distribution of transition radiation are dependent upon the dielectric constants of 
the two media and that transition radiation is not connected with changes in 
velocity of the charged particle. This behavior of transition radiation, particularly 
the angular distribution and polarization dependence on energy makes it a 
promising technology for use as a diagnostic for charged particle beams at low 
energy. 

Since its introduction, many investigations of the nature of transition 
radiation in such diverse environments as the boundary between two media, 
electromagnetic fields, plasmas, inhomogeneous media, and media in which the 
properties vary with time have been published, providing an exhaustive 
bibliography. Ginzburg and Tsytovich [Ref. 2], Ter— Mikaelian [Ref. 3], and 
Garibian [Ref. 4] provide brief overviews of some of the established concepts 
resulting from those investigations. Unfortunately, the calculations on the subject 
were always awkward in nature, since the formulas developed were very 
complicated. Many formulas developed independently of each other differed 
considerably from one publication to the next. The situation was exacerbated by 
the lack of specifics or standardization with respect to definitions of key concepts 
and adoption of universal conventions. Wartski [Ref. 5] attempted to rectify the 
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problem by returning to fundamental theories, building upon basic hypotheses, and 
emphasizing the conventions adopted by early pioneers in this field. Rule and 
Fiorito et al. [Ref. 6—8] refined and extended Wartski's work to develop analytical 
solutions concerning the components of radiation intensity. Much of the analysis in 
this thesis is based upon those analytical solutions. 

It is necessary to first discuss the physical processes involved in which 
transition radiation occurs. The simplest case to consider is that of a charged 
particle moving at a uniform velocity. which encounters a boundary between two 
media. The media are characterized by their respective dielectric permittivities, 
and e 9 . The dielectric permittivities are complex functions of the refractive index 
n and the absorption coefficient k of the medium. The motion of the charged 
particle creates, by definition, an electrical current which in turn generates 
associated electric and magnetic fields. Since the particle is assumed to travel in a 
uniform motion in which it does not experience acceleration, it does not emit 
radiation until encountering the boundary between the media. 

At least three types of radiation are emitted as the particle transits through 
the interface between the media. These are Cherenkov, bremsstrahlung, and 
transition radiation. Transition radiation is closely related to Cherenkov radiation 
but exhibits different properties. Bremsstrahlung radiation is an entirely separate 
mechanism involving particle acceleration. The form of transition radiation may be 
found by examining the fields involved. The electric and magnetic field in each 
medium are described by the solutions to the Maxwell equations for a point charge 
moving at constant velocity. However, the electric and magnetic field components 
in each media fail to satisfy the continuity requirements across the interface. 
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The solutions of homogeneous Maxwell equations must be added to those of the 
electric and magnetic fields in each medium to meet the boundary conditions. The 
required radiation fields described by these solutions meeting the boundary 
conditions are the transition radiation. [Ref. 4,9] 

Obviously, finding the solution for the transition radiation fields is quite 
complicated. Ter— Mikaelian [Ref. 3] and Wartski [Ref. 5] provide an excellent 
account of the full derivation. The results of those derivations will be presented and 
analyzed in this thesis so as to promote a better understanding of transition 
radiation and its potential applications for beam diagnostics. 

As stated earlier, the simplest case is that of a charged particle moving at 
uniform velocity encountering a boundary between media. For purposes of this 
Study, a vacuum to medium transition was assumed in which 1 and can be 
written simply as c. The medium encountered is a thin metallic foil. Transition 
radiation may be observed from the front face of the foil and is called backward 
transition radiation since it is reflected back from the foil. Transition radiation is 
also emitted in the forward direction from the back side of the foil and is thus 
termed forward transition radiation. Observation of fonvard transition radiation is 
difficult because of its proximity to the beam of charged particles. The backward 
transition radiation, however, may be diverted from the beam by reflecting it at an 
angle. For this reason, the beam of charged particles may be made to encounter the 
metallic foil at an oblique angle of 45 degrees. This has historically been the 
preferred angle of incidence because, since the angle of reflection equals the angle of 
incidence, it allows the observation equipment to be set up 90 degrees from the 
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beam axis along the axis of reflection. This configuration permits the study of 
backward transition radiation apart from the beam while keeping calculations 
relatively simple. 

The mean angle at which optical transition radiation appears is the angle of 
specular reflection. If the beam has an angular divergence then the charged 
particles have different angles of incidence upon the transition radiation foil. If the 
angle of a particular particle is a with respect to the beam axis, then the optical 
transition radiation will appear at an angle a to the axis of specular reflection. 
Figure 1 demonstrates the definition of the particle angle a with respect to the 
beam axis, angle of specular reflection T, and angle of observation 0. 




Figure 1. Definition of the particle angle a with respect to the beam axis, angle of 
specular reflection T, and the observation angle 0 with respect to the angle of spec 
reflection for a particle encountering a boundary at an oblique angle, [lief. 5,9]. 
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The interest here is the intensity of the radiation at any point in the plane of 
observation as a function of the observation angle 6. The angular information of 
transition radiation obtained from this is valuable as a potential diagnostic for 
charged particle beams. The TR3 program used in this study was based upon the 
analytical solution of this single foil case developed by Rule and Fiorito [Ref. 6-8] 
for this purpose. This program extends Rule and Fiorito's solution into three 
dimensions, thus using the advantage of three-dimensional graphics to aid in 
interpreting the results. 

Recall that the motion of the charged particle induces electromagnetic fields in 
the medium and that radiation is emitted when the particle encounters a boundary 
between two media of dissimilar dielectric constants. When that transition from 
one medium to another occurs, the electromagnetic fields induced by the particle 
motion as well as the fields induced by the transition must adjust to the change in 
dielectric permittivity. The distance over which this adjustment occurs is known as 
the coherence length. (It is also sometimes referred to as the formation length). 
The definition of this length is a bit arbitrary, but is usually chosen so that the 
radiation emitted at the interface adds coherently with the fields generated by the 
particle. It has been found that a phase difference of one radian is the condition for 
the rays to become coherent. The fields at much greater distances may be 
considered to be pure radiation fields. 

Coherence length is related to both the particle velocity and to the phase 
difference between the particle's fields and the radiation fields from the transition. 
The coherence length in the medium can be expressed [Ref. 5] as 
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where /? is the ratio of the particle velocity to the speed of light, u is the angular 
frequency ( =2ttc/ A), 6 is the angle of observation, and A is the wavelength of the 
radiation in the medium. The coherence length in a vacuum may be found by 
setting v^=l- 

Coherence length takes on a special significance when considering a particle 
encountering more than one boundary as with the two— foil interferometer developed 
by Wartski [Ref. 5]. This interferometer consists of two parallel foils set at an angle 
of 45 degrees to the trajectory of the charged particle beam. Particles traversing the 
foils emit radiation from the back of the first foil that is reflected off the front 
surface of the second foil. Transition radiation is also produced by the transit 
through the second foil. The radiation emitted from the front of the second foil 
forms an interference pattern with the radiation from the back of the first foil 
(Figure 2). The interference pattern is centered around the axis of specular 
reflection. The phase difference between the transition radiation emitted from the 
first and second foils is 

<J) = (2ttL/A/?)( 1— /?cos 6) = L/L y , (1.2) 

where L is the separation between the foils and L y is the formation length in a 
vacuum. This phase difference should be an integer multiple of one radian for the 
transition radiation from the two foils to add coherently. That is, the separation 
between the foils should be an integer multiple of the coherence length in a vacuum. 
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Figure 2. Diagram of t he Wartski's two— foil interferometer. [Ref. 5,9] 

Rule and Fiorito Ref 0—8] developed an analytical solution to the transition 
radiation emitted from Wartski's two— foil interferometer [Ref. 5]. The C011ER3 
program used in this study extended that solution into three dimensions for use with 
three-dimensional graphics to aid the analysis. 

Most current theories on transition radiation were built on the simplification 
of assuming a point value for both energy and frequency in order to evaluate 
complex integrals. In reality, however, energy and frequency extend over a finite 
range of values. The assumption may be responsible for some slight differences 
between theoretical and actual data. The capability of examining the effects of a 



t 



finite bandwidh of wavelength (and thus frequency) was built into the COHER3 
program. This thesis provides one of the few investigations on finite frequency 
bandwidth effects. 

The TR3 and COHER3 programs created an output compatible with the 
Three-dimensional graphics program, SON OF SURF3D. This graphics program 
was written by Don Gilbert of Dogstar Software as an extension of SON OF SURF. 
At the time of this writing, SON OF SURF3D was still under development and not 
yet ready for public release. This thesis was one of the first thorough testings of the 
SON OF SURF3D graphics program. As such, it may play a key role in 
determining whether the SON OF SURF3D program will be fully developed for 
commercial use. Use or discussion in this thesis of any comercially availabe product 
does not constitute endorsement. 

The purpose of this thesis was to promote a better understanding of transition 
radiation based upon the theoretical developments presented in references 1—9. It 
provides an investigation of the effects of finite bandwidth on radiation intensity 
distribution. The programs and graphics package used in this thesis could become 
powerful research tools in the development of TR as a diagnostic for charged 
particle beams. The behavior of the polarization and angular distribution of 
transition intensity in response to changes in specific parameters demonstrated in 
this thesis should provide a guide for further development of transition radiation 
theory for this purpose. By building upon basic concepts and examining the nature 
of the complex formulas involved, with the aid of three-dimensional graphics, this 
thesis proposes to clarify the concepts of transition radiation so that even the 
laymen could obtain a basic understanding. 
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Chapter II descibes the analytical solutions of Rule and Fiorito for transition 
radiation from the single foil and two— foil interferometer cases [Ref. 6—9] and the 
TR3 and COHER3 programs used in this study which extended those solutions to 
three dimensions. The chapter closes with a brief description of the three 
dimensional graphics program, SON OF SURF3D, used in conjunction with these 
programs. 

Chapter III contains the analysis of the effects of energy, beam divergence, and 
dielectric properties on the distribution of transition radiation intensity for the 
single foil case. An investigation into the effects of coherence length and optical 
frequency bandwidth as well as energy, beam divergence, and dielectric properties 
for the two— foil interferometer case forms the content of Chapter IV. 

The conclusions resulting from the analysis as well as recomendations for 
further study, hardware, and software are discusssed in Chapter V. For 

convenience, a users manual for the single foil and two— foil interferometer programs 
and for the three-dimensional graphics program was developed and included as 
Appendix A. Appendix 13 contains the listing for the single foil program TR3 and 
the two— foil interferometer program COHER3. 
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II. CALCULATIONAL APPROACH 



The core of the analysis of this study are analytical solutions developed by D. 
W. Rule et al. of the Naval Surface Warfare Center [Ref 6—9] for single foil and for 
two— foil interferometer transition radiation. A series of programs developed by D. 
W. Rule et al. around those analytical solutions were used in the study of transition 
radiation characteristics to produce theoretical two-dimensional profiles of relative 
intensity [Ref. 7—9]. 

Rule's programs were modified to extend the analysis to three dimensions and 
to add flexibility in setting data parameters during run time. The resulting 
program for the single foil case is called TR3 and an interferometer program which 
took into account coherence length and frequency bandwidth is called COHER3. 
Both programs were tailored specifically to be compatible with the 
three-dimensional graphics program, SON OF SURF3D, written by Don Gilbert of 
Dogstar Software [Ref. 10]. SON OF SURF3D was in the developmental stage, not 
yet ready for public release as an improved version of SON OF SURF, extending the 
capabilities of that program to plot three-dimensional surfaces from data files. A 
description of SON OF SURF3D follows the discussion of the theoretical basis of the 
TR3 and COHER3 programs. 

A. SINGLE FOIL 

The program TR3 was a modification of TRADS1, which was one of a series of 
single— foil transition radiation programs written by D. W. Rule [Ref. 8]. The 
TRADSl program calculated a two-dimensional intensity distribution as a function 
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of angle measured from the angle of specular reflection. The angle of specular 
reflection was chosen to be forty— five degrees as measured from the target plane and 
ninety degrees from the beam axis (Figure 1) for the reasons discussed in chapter 
one. The program assumes a Gaussian distribution of beam divergence angle and a 
vacuum to medium transition. The medium is characterized by a dielectric 
constant which is in general complex. All calculations were based on energy and the 
incremental values of the observation angle (6), measured from the angle of 
specular reflection. The TR3 program extended the analysis of the TRADS1 series 
of programs to three-dimensions in conjunction with the three-dimensional 
graphics program SOX OF SURF3D. 

The values of the parallel and perpendicular components of intensity as well as 
the total intensity were based on an analytical approach developed by D. W. Rule 
and associates [Ref. 7,8] in accordance with the theoretical treatment of transition 
radiation by Wartski [Ref. 5]. Using the small angle approximation of cos 0 and sin 
0. Rule [Ref. 7] showed that, for small angles, the intensity per unit frequency and 
solid angle could be expressed 



where ft is the velocity of the particle in units of the speed of light, Pjs the 
component of the velocity that is perpendicular to the observation plane, 7 is the 




( 2 . 1 ) 
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dimensionless Lorentz factor for energy, and ry and r^ are respectively the parallel 
and perpendicular Fresnel reflection coefficients. The Lorentz factor 7 and the 
velocity /? are defined 



7 = 



E 

9 

m o c 



and 





P = 




7“ 



v 

c ■ 



(2.3) 



(2.4) 



The parallel and perpendicular reflection coefficients ry and r^ are calculated 
in terms of the the observation angle 0 and dielectric constants of the media. As 
discussed earlier, a vacuum to medium transition was asssumed in which the 
dielectric constant e j of the vacuum is equal to one, and the dielectric constant of 
the medium may be written as e. Then rn and r are given by 



r 



ecos# — 



ccos 0 + 



• 2 a 

r-sin 0 



1 e— sin “0 
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( 2 . 6 ) 



Equations ( 2 . 1 ) and ( 2 . 2 ) contain only the two highest terms in the Lorentz 
factor 7 [Ref. 8]. Note that the odd— powered factor of 6 in In gives rise to an 
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asymetric pattern as the angle changes sign. The TR3 program preserves this 
asvmetric behavior by reversing the sign of the angles represented in the first 32 
columns in the 64 by 64 array used to calculate and organize the data. 

Rule's analytical calculations took into account the effect of beam divergence 
on the transition radiation intensity pattern. This was done by folding Gaussian 
distributions for the projected angles into equations (2.1) and (2.2), resulting in 



i„=^(2*r 3 /V 



W(Z)- 



7TCT 



X 



{|r||| 2 Re[(7+V?Za x ) 

1 1/2 

]— 2Re( r|| )ImW(Z)}, 



and 



T e 2 ^2/ 0 x— 3/2 -12 2, ,2 

\ i = r 0-(2,) '^yrlrj 

x Re[(r-4Z/CT v )W(Z)+ur 1 (2/i) 1 / 2 ] , 



(2.7) 



( 2 . 8 ) 



where Z is defined as (7 _1 +i0)/,/5<7 , a as the rms beam angle of divergence 
projected into the observation plane, and a as the rms angle of beam divergence 
perpendicular to this plane and containing the beam axis. Rule ignored the second 
term of equation (2.2) for equation (2.8) and it was also assumed that the rms beam 
divergence was small, i.e. o« q -1 . The function W(Z) is given in terms of the 
error function cj)(Z) in the non-standard form 



W(Z) = [l-4>(Z)]exp(Z 2 ). (2.9) 



The whole of the single foil program, TR3, was built around the analytic 
solutions embodied by equations (2. 7-2. 9). Total intensity was defined simply as 
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the sum of the perpendicular and parallel components. All intensities were 
calculated in terms of charge squared divided by the speed of light per unit 
frequency and solid angle and were normalized for plotting purposes. Actual values 
of intensities as defined here were stored in a file labelled RADOUT.DAT. Note 
that if the beam divergence is symmetric such that <7^= o then the total 
divergence angle is fio. This case was built into the single foil program. 

Rule's analytic solutions were extended into three dimensions by calculating 
the parallel and perpendicular components of radiation intensity and the total 
radiation intensity over an angular distribution centered around the angle of 
specular reflection. This angular range was broken down into a 64 by 64 array in 
which the intensities were calculated for each point in the array. The parallel 
component, perpendicular component, or total intensity for each point may be 
stored in the graphics output file for plotting. Each point in the array represents an 
angle and position measured from the axis of specular reflection. The axis was set 
at the center of the array in the thirty— second column of the thirty— second row of 
the array. In terms of data storage, the rows of the array correspond to records and 
the columns correspond to field elements. Each of the 64 records contain 64 field 
elements. Since the program takes several minutes to calculate the data, progress is 
indicated by a record number. 

B. INTERFEROMETER 

The interferometer program, COHER3, was based on the TEM series of 
transition radiation programs developed by D.W. Rule et al. that culminated in the 
two-dimensional program COHER1. The COHERl program was dveloped around 
Rule's analytical solutions for transition radiation in the case of Wartski's two— foil 
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interferometer. A finite optical frequencycy bandwidth was also built in, allowing 
for a study of the affects of finite bandwidth on transition radiation intensity 
distribution. The COHER3 program simply extended COHER1 to three 
dimensions, modifying it so as to allow flexibility in altering data parameters during 
run time, and tailoring the output to be compatible with the three-dimensional 
graphics program, SON OF SURF3D 

At the time of this writing, Rule and Fiorito had not published a full 
discussion of their analytical solution for the two— foil interferometer case. However, 
a discussion of the development of the analytical code for parallel component of 
intensity was published [Ref. 7] and formed the basis of the presentation here. The 
perpendicular component follows a. similar argument in development. 

Rule and Fiorito's analysis began with Wartski's [Ref. 5] development of the 
two— foil interferometer. Wartski showed that the parallel component of TR 
intensity per unit frequency and solid angle in this case could be written 



where u represents frequency, fl represents a solid angle, ip represents the angle of 
specular reflection, 0 represents the observation angle as measured from ip, and 
F (ip,0, J) represents a Fresnel coefficient which for a vacuum to medium transition is 
defined 




4jt^c ( 1 —0 c o s 6 ) ^ 




( 2 . 10 ) 



F(ip, 6.-') = i f ~ si — tcos {t±&} 



[f-sin ^(ip+6)]^^+ ecos(ip+9) 



( 2 . 11 ) 
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The exponential term is the interference term for the transition radiation 
patterns differing in phase by <j>. This phase difference is defined in equation (1.2) 
as the ratio of the interfoil spacing to the coherence length. For relativistic 
particles, with a Gaussian distribution of beam angles folded in and averaged over 
particle angles with respect to the beam a, the intensity per unit frequency and 
solid angle may be written 



The basis of the COIIER3 program is the analytical code developed by Rule 
and Fiorito to solve both equation (2.12) and the integral for the perpendicular 
component of intensity as well. The result for each component may be written 





( 2 . 12 ) 
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(2.13) 




X 




(2.14) 
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C. THE GRAPHICS PROGRAM 

The three-dimensional graphics program used in conjuction with the TR3 and 
COHER3 programs for the analytical study of transition radiation was SON OF 
SURF3D [Ref. 10]. This program was a derivative of a program written by Don 
Gilmore of Dogstar Software to explore the graphics capabilities of TurboHalo 3.0 
[Ref. 11]. SON OF SURF3D is able to plot three-dimensional surfaces from data 
files containing three-dimensional point information of up to 100 by 100 arrays. 

The data points stored in the output file for the graphics program SON OF 
SURF3D are automatically standardized relative to the axes according to the 
following formula [Ref. 10] 

Z' = ^ ^min ^cen^ . (2.15) 

TZ -Z . ) 
v max min' 

The TR3 and COHER3 programs store normalized real values of the parallel 
component, perpendicular component, or total intensity in an output file named by 
the user. For purposes of the analysis, the z— maximum value was set at 1.0 and 
both the z— minimum and z— center values were set at 0.0. Doing so sets the axis 
markers at the bottom of the data set in the three-dimensional plots and shows the 
intensties relative to a maximum intensity. An option was provided to normalize 
all intensities to the maximum intensity at 100 MeV. 
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Ill SINGLE FOIL TRANSITION RADIATION 



The TR3 program was run many times while changing one parameter each 
time in order to study the effect of that parameter on the distribution of 
transition radiation intensity. The effects of energy, beam divergence, 
scattering foil, and dielectric constant foil were studied in this manner. 

A. ENERGY 

The single foil case for transition radiation was calculated for energies of 1 
GeV, 500 MeV, and for 10 MeV to 100 MeV in increments of 10 MeV. The total 
intensity distribution as a function of angle was calculated over an angular 
range of —.05 to + .05 radians with respect to the angle of specular reflection. 
Several observations were made from the resulting data (Tables I, II) and 
associated three-dimensional plots (Figures 3—7). The data for Figures 3—7 
were normalized to the maximum value of intensity for the energy plotted. 
Therefore the maximum relative value is one (1.0), which is represented by the top 
of the z— axis. The units on the z— axis are thus the dimensionless ratio of intensity 
to maximum intensity. The angular distribution was over 0.05 radians as measured 
from the z— axis in either direction along the x and y axes. The endpoints on the x 
and y axes coincide with .0625 radians. 

The most obvious observation was that the shape of the plots varied 
greatly from one energy to the next when using the fixed parameters described 
above. However, the plots appear nearly identical when plotting the data over a 
relative scale such as a multiple of the predicted angle of maximum intensty. 
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Such is the case for Figures S— 10 in which the data was calculated over an angular 
distribution of five times the predicted angle of maximum intensity for the energy 
plotted, i.e. | 0\ < 5/7. 

At low energy the plot was a gently curved cone shape centered around 
the z— axis with the opening extending almost the full angular width of the 
plot. The angular separation measured from the angle of specular reflection 
(z— axis) was defined as theta (Figure 1 of chapter I). The angle between the 
z— axis and the sides of the cone form a particular value of theta. As the 
energy was increased, the cone became more sharply defined and the angle 
steadily decreased. Wartski [Ref. 5] showed this phenomenon was predicted in 
theory from the equation for transition radiation intensity per unit frequency 
and solid angle in this case for a vacuum to medium transition. 



( 3 . 1 ) 

where q is the charge on the particle, c is the dielectric constant of the medium, 
and the fresnel coefficients r.. and fn are defined as 




2 



r 



II 




( 3 . 2 ) 
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Figure 3. Single foil transition radiation at 10 Mev. Note the overall cone shape of 
the plot. The sides of the cone describe the angle at which the maximum intensity 
occurs relative to the axis of specular reflection (Z-axis). The X-Y plane depicts 
the angular distribution in radians. The data extends to from -.05 to +.05 radians 
on each axis. The endpoints on the X and Y axes correspond to ±.0625 radians, 
fhe units measured along the Z— axis are the dimensionless ratio of intensity to 
maximum intensity. The top of the Z— axis corresponds to the maximum value of 
the data. 
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Figure 4. Single foil transition radiation at 40 Mev. The cone shape still exists in 
the center of the plot, but the angle describing the position of maximum intensity is 
considerably smaller. The X— Y plane depicts the angular distribution in radians. 
The data extends to from —.05 to +.05 radians on each axis. The endpoints on the 
X and Y axes correspond to ±.0625 radians. The units measured along the Z— axis 
are the dimensionless ratio of intensity to maximum intensity. The top of the 
Z— axis corresponds to the maximum value of the data. 
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Figure 5. Single foil transition radiation at 70 Mev. The original shape of a cone is 
nearly lost in the center of the plot. The angle describing the position of maximum 
intensity is again much smaller. The X-Y plane depicts the angular distribution in 
radians. The data extends to from -.05 to +.05 radians on each axis. The 
endpoints on the X and Y axes correspond to ±.0625 radians. The units measured 
along the Z-axis are the dimensionless ratio of intensity to maximum intensity. 
The top of the Z— axis corresponds to the maximum value of the data. 
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Figure 6. Single foil transition radiation at 100 MeV. The only remaining evidence 
of the original cone appears to be the rim occurring at the top of the plot around the 
Z-axis. The apparent asymmetry of the rim is due to the plotting algorithm and is 
not a physical reality. The X-Y plane depicts the angular distribution in radians. 
The data extends to" from -.05 to +.05 radians on each axis. The endpoints on the 
X and Y axes correspond to ±.0625 radians. The units measured along the Z— axis 
are the dimensionless ratio of intensity to maximum intensity. The top of the 
Z— axis corresponds to the maximum value of the data. 
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Figure 7. Single foil transition radiation at 1 GeV. The data appears to converge to 
a point below the top of the Z— axis which should represent the maximum value oof 
the data. In fact, the maximum value occured on the axis of specular reflection 
represented by the Z— axis. Furthermore, maximum value of the data occurred as a 
point which is at the top of the Z-axis. The X-Y plane depicts the angular 
distribution in radians. f l he data extends to from —.05 to +.05 radians on each axis. 
The endpoints on the X and Y axes correspond to ±.0625 radians. The units 
measured along the Z— axis are the dimensionless ratio of intensity to maximum 
intensity. T he top of the Z— axis corresponds to the maximum value of the data. 
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Figure 8 . Single foil transition radiation at 30 Mev. Compare with the plot for 
single foil transition radiation at 40 MeV (Figure 4). In this case data extends to 
five times the predicted value of maximum intensity, 1/7 in either direction on the 
X and Y— axes. The Z— axis still measures the relative intensity of the data with the 
top of the Z-axis corresponding to the maximum intensity. 
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Figure 9. Single foil transition radiation at 60 Mev. The plot is nearly identical to 
the plot for 30 MeV (Figure 8), which is plotted on the same relative scale. That is, 
the data extends to five times the predicted value of maximum intensity, 1/7 
either direction on the X and Y-axes. The Z-axis measures the relative intensity ol 
the data with the top of the Z-axis corresponding to the maximum intensity. 
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Figure 10 . Single foil transition radiation at 90 Mev. Compare with the plot for 30 
and 60 MeV (Figures 8,9), measured on the same relative scale. The plots appear to 
be identical. The shape of the single foil transition radiation pattern is independent 
to the energy. Intensity increases with energy increases while the angle at which the 
maximum intensity occurs decreases. The data extends to five times the predicted 
value of maximum intensity, 1/7 in either direction on the X and Y— axes. The 
Z— axis still measures the relative intensity of the data with the top of the Z— axis 
corresponding to the maximum intensity. 



27 



f„ = 



2 e cos 9 



fcos 9 4- /e-sin 2 $ 



The fresnel coefficients must satisfy the continuity relationship: 

1+r n - f n 

For metallic media in the visible region, | c | >>1 , and equation 
(3.1) can be rewritten 



1 = 



0 2 • 2 n 
1 e sin 9 



4tt 2 c (1-/?cos0) 2 



Since /? ^ 1 — ( 1 / 2 7 ) for the extremely relativistic case 



I = 



2 e^ 



sin 2 # 



4x 2 c [1 -cos#+(1/2 7 2 )cos#] 2 



A second approximation was made by considering the angl 

2 2 

sufficiently small such that sin 9 and cos 9 could be replaced b; 
9 

(1 -(T / 2 ), respectively, leaving 
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(3.5) 
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Again, the preceeding equations were for the specific case at hand of a 



vacuum to medium transition. The factor |i|| , the reflectance of the 
medium, can be ignored for a medium to vacuum transition. Equation (3.7) 
showed the relationship of radiation intensity to energy and the angle 
theta. The nature of this relationship was shown by differentiating 
equation (3.7) with respect to theta: 



dl 



d 6 ~~c 



2 [2»(l/ 7 2 - e 2 )} 

(£ + 1 / -) 2 ) 3 



(3.8) 



This equation demonstrated that the extrema occur for the theta equal 
zero direction and for theta equal 1/7. Using these values for theta in equation 
(3.7) showed that the radiation intensity is zero in the theta equal zero direction 
and is at a maximum for the theta equal 1 / 7 direction. This angle of maximum 
radiation intensity decreased as energy increased. A look at Figures 3—7 
showed that the relationship was qualitatively true. A close look at the 
data (Table I) . in which the angle measured where maximum intensity actually 
occurred was compared with 1/7 demonstrated excellent correlation. 

The energy, the Lorentz factor 7, predicted angle of maximum intensity 
1/7. and measured angle of maximum intensity 0 ^, for a beam divergence of .001 
radians were tabulated in Table I. The maximum value of the total intensity 
ITOT, maximum value of parallel intensity IPAR, maximum value of 
perpendicular intensity IPERP, the value of parallel intensity at maximum total 
intensity IPAR(# m ), and the value of perpendicular intensity at maximum total 
intensity IPERP(# m ). were tabulated in Table II for comparison. Intensity 
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was calculated per unit frequency and solid angle in terms of the square of 
the charge divided by the speed of light. Energy was listed in units of MeV, 
angles were listed in units of radians. Additionally, parallel, perpendicular, 
and total intensities at the angle of maximum intensity as well as the maximum 
of the perpendicular intensity versus energy were plotted in Figure 11. The 
plot revealed that total and parallel intensity increased as a power 
function of energy. 



TABLE 1. ACTUAL VS. ANGLE OF MAXIMUM INTENSITY 



Energy 

(MeV) 


7 


1/7 

(radians) 


Actual 0 
(radians) 


10 


20.57 


.0486 


.0487 


20 


40.14 


.0249 


.0250 


30 


59.71 


.0168 


.0169 


40 


79.28 


.0126 


.0126 


50 


98.85 


.0101 


.0102 


60 


118.4 


.0084 


.0087 


70 


137.0 


.0073 


.0072 


SO 


157.6 


.0064 


.0065 


90 


177.1 


.0057 


.0058 


100 


196.7 


.0051 


.0051 


500 


979.5 


.0010 


.0000 



It was interesting to note that the maximum intensity at .5 GeV and 1.0 
GeV occurred at the theta equal zero direction. A look at the data in Table II 
provides a clue as to why this occurred. As energy increased, the perpendicular 
component of intensity provided a greater proportion of the total intensity. It was 
found that the perpendicular component of intensity began to dominate over the 
parallel component at energies above 200 MeV. At .5 GeV and 1.0 GeV, the 
perpendicular component of intensity dominated over the parallel component of 
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Figure 11. A plot of the components of intensity versus energy. It is apparent that 
the perpendulai component at the angle where the maximum intensity occuis is 
negligible at lower energies. However, the maximum value of the perpendicular 
component can not be ignored. This maximum value occurs on the axis of specular 
reflection, and thus tends to fill up the cone shape seen in Figures 3-5. 
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TABLE II. INTENSITIES 



ENERGY 

(MeV) 



ITOT IPAR IPERP IPAR(0 m ) IPERP ( 0 m ) 



2, , 
e / c 




r 2, , 

e /c 




2/ , 
e / c 




r 2, i 
e /c 




r 2/ , 

e / c 


doxR2_ 




dodfi_ 




doxlfl. 








dodO_ 



10 


9.571 


20 


36.51 


30 


80.83 


40 


142.5 


50 


221.7 


60 


318.1 


70 


432.3 


80 


564.1 


90 


713.5 


100 


881.4 


500 


5S230 


1000 


467300 



9.567 


0.171 


36.45 


.2479 


80.52 


1.209 


141.6 


3.740 


219.4 


8.979 


313.5 


18.34 


423.2 


33.53 


549.3 


56.42 


689.9 


89.13 


843.0 


133.9 


14230 


44740 


52090 


415200 



9.567 


.0043 


36.45 


.0618 


80.52 


.3036 


141.6 


.9617 


219.4 


2.296 


313.5 


4.580 


423.2 


9.106 


549.1 


15.01 


689.9 


23.63 


843.0 


38.36 


13480 


44740 


52090 


415200 



intensity in contributing to the total intensity. The components of radiation 



intensity behave quite differently from each other and contribute distinct profiles to 
the total intensity distribution. Wartski [Ref. 5], (Figure 12) showed that the 
perpendicular and parallel components of radiation intensity could be 



written: 
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47r 2 c 
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(3.10) 



where n, n', n", are as depicted in Figure 12 and 0 ^ and 0 ^ are the 
projections of the r/c vector on the plane of observation. Wartski [Ref. 5] 
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z 

* 




Figure 12. The the unit n, n', ii", and $ vectors in the plane of observation as 
depicted in Wartski [Ref. 5]. 



showed that for ultra— relativistic particles, where it is assumed that y/7» 1 
in a vacuum to medium transition, the second term in equations (3.9) and 
(3.10) dominates as /)-ii becomes close to one. In this case 
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(3.12) 
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The maximum value of 1^ does not occur at the same position as the 
maximum of total intensity except at very high energies where the perpendicular 
component of intensity becomes dominate and thus determines the location of the 
maximum total intensity. It was found that the energy where this occurs is in the 
neighborhood of about 200 MeV. A look at the perpendicular reflection coefficient 
reveals that the maximum of the perpendicular component of radiation intensity 
occurs when theta (0),the angular separation from the observation angle (z— axis) is 
zero 



r 

-L 



/ o — 

cos 0 — V e — s in “ 0 

cos 0 + V e — s i n ^ 6 



1 



- Vc 



1 + € 



0=0 



(3.13) 



For non— relativistic particles, the perpendicular component of the 
radiation intensity could be shown to be essentially zero while the parallel 
component of radiation intensity could be written 



t e q2 • 2 n 

’ir^z 511 ' K 
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1 + r. 



41 



42 9 

= e /r 

2„ 

7T C 



6-1 



sin^ 0 cos ^ 0 
ccos 0 +Jc-sin'0 I 2 



(3.14) 



The 3 in equation (3.14) indicates that the radiation produced by a 

Zj 

non— relativistic particle for an oblique angle is the same as that for 
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normal incidence at a velocity of ft rather than ft. It was clear from 
equation (3.14) that the parallel component of radiation intensity for a 
non— relativistic particle is proportional to the square of the energy as is 
the total intensity. 

In the ultra— relativistic energy region, Iy remains proportional to 
the square of the energy. Equation (3.12) indicates that I is also proportional to 
the square of the energy for lower energies. The contribution of I to the total 
energy is negligible at lower energies but becomes quite significant as energy 
increases, and increases asymptotically as /3-n approaches one. Figures 13—15 depict 
the perpendicular component of radiation intensity at 100, 250, and 500 MeV. 

Analysis of equation (3.7) revealed that total intensity was 
proportional to the square of the energy. Figures 16—18 dramatically showed 
the dependence of radiation intensity on energy. The intensity 
distribution for these figures was scaled to the maximum intensity value of 
100 MeV. The energies for these figures were 50 MeV to 90 MeV in increments 
of 20 MeV. 

Finally, a closer look at Figures 8-10 reveals that the plots aren't quite 
identical in appearance. In each case the intensity distribution was calculated over 
an angular distribution of five times the predicted angle of maximum intensity of 
l/'z for the energies of 30. 60. and 90 MeV. Closer inspection revealed that the apex 
of the cone decreased in depth as energy increased. Obviously, the apex of the cone 
was not zero intensity at higher energies. This was due to the contribution of the 
perpendicular component of intensity to the total intensity distribution. 
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Figure 13. The perpendicular component of intensity at 100 MeV. The data 
extends to five times the predicted value of maximum intensity, 1/7 in either 
direction on the X and Y— axes. The Z— axis measures the relative intensity of the 
data with the top of the Z— axis corresponding to the maximum total intensity. 
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Figure 15. The perpendicular component of intensity at 500 MeV. The data 
extends to five times the predicted value of maximum intensity, 1/7 in either 
direction on the X and Y-axes. The Z-aocis measures the relative intensity of the 
data with the top of the Z— axis corresponding to the maximum total intensity. 
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Figure 16. Single foil transition radiation intensity at 50 MeV. The data extends to 
five times the predicted value of maximum intensity, 1/7 in either direction on the 
X and Y-axes. The Z-axis measures the intensity of the data, relative to the 
maximum intensity obtained for 100 MeV using the default values of the program. 
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Figure 17. Single foil transition radiation intensity at 70 MeV. The data extends to 
five times the predicted value of maximum intensity, 1/7 in either direction on the 
X and Y— axes. The Z— axis measures the intensity of the data relative to the 
maximum intensity obtained for 100 MeV using the default values of the program. 
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Figure 18. Single foil transition radiation intensity at 90 MeV. The data extends to 
five times the predicted value of maximum intensity, 1/7 in either direction on the 
X and Y-axes. The Z-axis measures the intensity of the data relative to the 
maximum intensity obtained for 100 MeV using the default values of the program. 
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TABLE III, MULTIPLES OF 1/7 



Energy 

(MeV) 



7 



1 h 5 /7 



(radians) (radians) (radians) 



10 

30 

60 

90 



20.57 

59.71 

118.4 

177.1 



.0486 

.0168 

.0084 

.0057 



.2431 

.0837 

.0422 

.0282 



.0487 

.0168 

.0087 

.0058 



In summary, the angle of maximum intensity, total intensity, perpendicular 
componenet of intensity, and the parallel component of intensity were found to be 
very sensitive to changes in energy. The perpendicular component of intensity was 
found to be rather insignificant at energies below 50 MeV but quite dominate at 
energies above 200 MeV. The angle of maximum intensity, total intensity, and 
parallel component of radiation intensity would be the most useful parameters for 
diagnosic purposes of beam quality at lower energies. However the angular 
information is lost above 200 MeV when the perpendicular component of intensity 
becomes dominate. 

B. BEAM DIVERGENCE 

The effect of beam divergence upon the transition radiation intensity 
distribution may be studied by assuming a Gaussian distribution of divergence 
angles. The Gaussian distribution of angles may be approximated by taking the 
root mean square average value and applying it to the transition radiation intensity 
equations (2.7) and (2.8) discussed in Chapter II. 
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Recall that when writing the code for the single foil case, Rule [Ref. 7] 
programmed in the assumption that a - a = a such that the total divergence 

X v 

angle is 4a. In this case, equations (2.7) and (2.8) may be rearranged to bring out 
the dependence on beam divergence more clearly. 
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It is now clear that the parallel component of intensity should decrease with 
an increase in beam divergence. On the other hand, any increase in beam 
divergence should bring about a corresponding increase for the perpendicular 
component of intensity. Recall that the maximum intensity of the perpendicular 
component occurs on the axis of specular reflection. The maximum of the parallel 
component, however, occurs away from that axis at an angle equal to the inverse of 
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the Lorentz factor. Suppose the beam divergence would be permitted to become 
large enough for the perpendicular component to become comparable with, or even 
dominant over the parallel component. The angle at which the maximum intensity 
occurs would shift, eventually ocurring on the axis of specular reflection whenever 
the perpendicular component dominates over the parallel component. The angular 
information provided by the parallel component would then be lost. Both 
components of intensity increase as a power function of energy so beam divergence 
effects should be more pronounced at higher energies. 

The data in Table IV and in Figures 19—28 bear out the correctness of this 
analysis. Table IV contains a comparison of the total, parallel, and perpendicular 
intensities for various energies and rms beam divergence angles. Energy is listed in 
MeV, beam divergence, labelled Beam Div., is listed in radians. The intensities are 
per unit frequency and solid angle in terms of charge squared divided by the speed 
of light. The values for the parallel and perpendicular components of intensity are 
the maximum values obtained for that component. Since the maxima of the 
components of intensity do not occur in the same position, the sum of the values 
listed will not agree with the value of the total intensity. The angle 0 at which 
the maximum intensity occurs is in radians. It is included here to study the 
behavior of this parameter as the parallel and perpendicular components of intensity 
adjust roles to changes in beam divergence. Note that an increase in beam 
divergence is accompanied by a corresponding increase in the total and perpendicuar 
intensities and a decrease in the parallel component of intensity. Whenever the 
beam divergence was large enough for the perpendicular component to dominate 



44 



over the parallel component, the maximum intensity occurred on the axis. This can 

be seen in Table IV in that whenever the value of IPERP was larger than that of 

IPAR. the value of 9 was zero. 

m 

The data for Figures 19-28 was normalized to the maximum intensity for the 
energy plotted. The top of the z— axis corresponds to the value of the maximum 
intenstv. The horizontal width of the data in the plots of Figures 19-26 was fixed 
at 0.05 radians as measured in either direction along the axis from the z— axis. The 
ends of the X and Y axes correspond to a value of .0625 radians. The horizontal 
width of the data in the plots of Figures 27 and 28 was set at a value of five times 
the predicted angle of maximum intensity. This option was chosen to preserve 
clarity at high energies. 

The effects of beam divergence on the intensity distribution are clearly seen in 
Figures 19—28. Figures 19—23 depict the effects of increasing beam divergence at 10 
MeY. A general flattening of the intensity distribution seems to occur. However, 
recall that the predicted angle of maximum intensity for this energy is .049 radians 
while the data is calculated out to .05 radians. Thus Figures 19—23 provide a close 
look at the region where the perpendicular component of intensity becomes more 
and more dominant. The perpendicular component of intensity appears to become 
completely dominant when the beam divergence was increased to just over .0250 
radians. 

Figures 24-26 depict the shape of a gaussian distribution at an energy of 40 
MeY as the beam divergence is increased from 0.0010 to 0.0060 radians. In this case 
the predicted angle of maximum intensity at .0126 radians is much less than the 
.0500 radians over which the data is calculated. Therefore, the effects of beam 
divergence on the parallel, as well as the perpendicular, component of intensity can 
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Figure 19. Single foil transition radiation pattern at 10 MeV for an rms beam 
divergence angle of 1.000E— 3 radians. The data extends to .05 radians in the X— Y 
plane. The top of the Z-axis corresponds to the maximum intensity of the data. 
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Figure 20. Single foil transition radiation pattern at 10 MeV for an rms beam 
divergence angle of 1.000E-2 radians. The depth of the cone is much shallower than 
in the previous figure. The data extends to .05 radians in the X— \ plane. The top 
of the Z— axis corresponds to the maximum intensity of the data. 
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Figure 21. Single foil transition radiation pattern at 10 MeV for an rms beam 
divergence angle of 1.250E-2 radians. Note that the depth of the cone continues to 
decrease with an increase in beam divergence. The data extends to .05 radians in 
the X-Y plane. The top of the Z-axis corresponds to the maximum intensity of the 
data. 
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Figure 22. Single foil transition radiation pattern at 10 MeV for an rmS beam 
divergence angle of 2.000E — 2 radians. What s left of the cone appears to e 1 e 
more than a dimple. Note that the angle at which the maximum intensity occurs fo 
an energy of 10 MeV is .049 radians while the data presented here extends to .05 
radians in the X-Y plane. The region depicted then coincides with the region m 
which the perpendicular component of intensity becomes more and more dominant 
over the parallel component. The top of the Z-axis corresponds to the maximum 
intensity of the data. 
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Figure 23. Single foil transition radiation pattern at 10 MeV for an rms beam 
divergence angle of 2.500E— 2 radians. All the angular information of the TR 
pattern has been completely washed out by the effects of beam divergence. The 
data extends to .05 radians in the X-Y plane. The top of the Z-axis corresponds to 
the maximum intensity of the data. 
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Figure 24. Single foil transition radiation pattern at 40 MeV for an rms beam 
divergence angle of 1.000E— 3 radians. The maximum intensity occurs at an angle of 
.0120 radians while the data, presented extends to .05 radians in the X— Y plane. 
The top of the Z— axis corresponds to the maximum intensity of the data. 
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Figure 25. Single foil transition radiation pattern at 40 MeV for an rms beam 
divergence angle of 5.000E-3 radians. Note that the rim around the cone is more 
rounded thatn in the previous figure and that the depth of the cone has significantly 
decreased. The increase in the perpendicular component, which is at a. maximum on 
the axis, is responsible for the decreasing depth of the cone. The data extends to .05 
radians in the X-Y plane. The top of the Z-axis corresponds to the maximum 
intensity of the data. 



52 




Figure 26. Single foil transition radiation pattern at 40 MeV for an rms beam 
divergence angle of 6.000E— 3 radians. The angular information has been nearly 
washed out by the beam divergence. The perpendiular component of intensity is 
nearly equal to the parallel component of intensity. The data extends to .05 radians 
in the X— Y plane. The top of the Z— axis corresponds to the maximum intensity of 
the data. 
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Figure 27. Single foil transition radiation pattern at 100 MeV for an rms beam 
divergence angle of 1.500E-3 radians. The data extends to five times the predicted 
angle of maximum intensity, i.e. ±5/7 in the X— Y plane. This option was chosen to 
preserve clarity at this high energy. The top of the Z-axis corresponds to the 
maximum intensity of the data. 
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Figure 28. Single foil transition radiation pattern at 100 MeV for an rms be.am 
divergence angle of 3.000E— 3 radians. All the angular information has been washed 
out by the effects of beam divergence. Note the apparent gaussian distribution ol 
the data. The data extends to ± 5/7 radians in the X— Y plane. I he top ol the 
Z— axis corresponds to the maximum intensity of the data. 
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TABLE IV, BEAM DIVERGENCE VS INTENSITIES 



Energy 


Beam Div. 




ITOT 




IPAR 




IPERP 0 


(MeV) 


(radians) 




'e 2 / c' 




re 2 /cl 




re 2 /c] 


(radians) 




daxHl_ 




doxlfl 




daxID 


10 


1.000E— 3 


9.571 


9.567 


1.711 E- 


-2 .0487 


10 


2.000E— 3 


9.571 


9.554 


6.826E- 


-2 .0485 


10 


3.000E— 3 


9.573 


9.534 


0.1529 


.0488 


10 


4.000E-3 


9.574 


9.506 


0.2703 


.0488 


10 


5.000E— 3 


9.578 


9.471 


0.4193 


.0488 


10 


6.000E— 3 


9.582 


9.428 


0.5987 


.0490 


10 


7.000E— 3 


9.588 


9.378 


0.8069 


.0490 


10 


8.000E— 3 


9.597 


9.322 


1.042 


.0490 


10 


9.000E— 3 


9.609 


9.259 


1.304 


.0491 


10 


1.000E— 2 


9.625 


9.191 


1.589 


.0490 


10 


2.000E— 2 


10.37 


8.318 


5.425 


.0423 


10 


2.500E— 2 


12.00 


7.853 


7.781 


.0178 


10 


3.000E— 2 


15.11 


7.424 


10.29 


0.000 


40 


2.500E— 3 


143.2 


137.2 


22.07 


.0126 


40 


5.000E— 3 


152.5 


124.9 


75.96 


.0114 


40 


6.000E— 3 


166.4 


119.5 


102.4 


.0082 


40 


7.500E— 3 


214.8 


111.9 


144.9 


0.000 


100 


5.000E— 4 


877.8 


868.7 


35.21 


.0051 


100 


1.500E-3 


895.2 


810.2 


281.1 


.0050 


100 


2.000E— 3 


937.1 


770.0 


461.7 


.0046 


100 


3.000E— 3 


1310 


690.3 


8S1.6 


0.000 



be seen. In Figures 24—26 the parallel component of intensity increasingly 
dominates the distribution at and near the z— axis while the parallel component of 
intensity continues to dominate the intensity distribution farther out from the axis. 
The overall intensity distribution approximates a Gaussian distribution once the 
beam divergence is large enough for the perpendicular component to provide the 
maximum overall intensity on the axis of specular reflection. The point at which 
this occurs for 40 MeV is for an rms beam divergence angle of just over 0.060 
radians. For 100 MeV, this point occured between 0.0020 and 0.0030 radians. 
Figures 27 and 28 display the intensity distribution at 100 MeV for rms beam 
divergence angles of 0.0015 and 0.0030 radians. 
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In summary, increased beam divergence increases the perpendicular 
component of intensity while decreasing the parallel component. The angle of 
maximum intensity decreases as the parallel and perpendicular components adjust 
in response to the increase in beam divergence. This can result in the eventual loss 
of all angular information as the intensity distribution approximates a Gaussian 
distribution. The effect is more pronounced at higher energies and thus a beam 
divergence that may be negligible at low energies may have significant effect at 
higher energies. Beam divergence has a detrimental affect on beam quality, 
especially at higher energies and should be minimized to preserve angular 
information from transition radiation. 

C. DIELECTRIC CONSTANT 

As has been stated before, transition radiation occurs when a charged particle 
in motion encounters a boundary between media. The intensity of the radiation 
emitted is dependent upon the dielectric properties of the media involved. The 
majority of the study of transition radiation here has assumed a vacuum to medium 
transtion. That assumption simplifies the analysis and the formulas involved. Now, 
however, it is more useful to study the equation for the total intensity per unit 
frequency and solid angle in which neither medium is assumed to be a vacuum. 

The classic case, as developed by Frank and Ginsberg [Ref. l],.is to assume 
that the particle travelling through a medium which is characterized by adielectric 
constant transits an interface to another medium characterized by a dielectric 
constant t. r The trajectory of the particle is assumed to be normal to the 
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interface. Frank and Ginsberg [Ref. 1] developed expressions for the transition 
radiation intensity in both media. The intensity observed in the first medium per 
unit frequency and solid angle at an angle 0 ^ from the the normal is 



yf77 sin^.cos ^6, 2 

h = s — 

7 r c 

( 1 —/^ ) 

(1— (?“ t^cos^ 0^)(1 — (ij e 2 -e l s * n ^^l) ( e 9 C0S ^l'l'V^f^fo — 



(3.16) 



The intensity per unit frequency and solid angle observed in the second 
medium at an angle 0 9 to the normal may be found from the above expression by 
substituting — (3 for /? and interchanging the subscripts 1 and 2. In either case it is 
clear that the intensity observed in the medium is proportional to the square of the 
difference of the dielectric constants of the two media. 

Actually the dielectric constant isn't really a constant at all. It is really a 
slowly— varying function of wavelength that may be considered a constant when 
considering short ranges of wavelength. The dielectric constant is complex, 
consisting of both a real and imaginary part that are related to the refractive index 
n and the absorption coeffiecient k of the medium. Letting fl represent the real 
part of the dielectric constant e and e2 the imaginary part, the following relations 
define the dielectric constant; 



e = el + ie2 , 


(3.17) 


1 

Cd 

Sh 

II 

t— H 

S4J 


(3.18) 
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c2 = 2nn . 



(3,19) 



Values for the refractive index and absorption coefficient at discrete 
wavelengths for some metals were found in the Handbook of Chemistry and Physics 
[Ref. 12]. The available information on the dielectric constants of the metals 
limited the data analysis on this subject to three significant figures. TABLE V lists 
the optical wavelength, refractive index n, absorption coefficient k, real part of 
the dielectric constant el, and the imaginary part of the dielectric constant e2, for 
aluminum Al, gold Au, and silver Ag. Note that the value listed the real part of 
the dielectric constant is the absolute value of the result of equation (3.18). 

Some of the results obtained by this method required verification. The 
formulas listed in equations (3.18) and (3.19) were also found in Wartski [Ref. 5]. 
and in Bennett and Bennett (Ref. 13] which also provided a table of optical 
constants at infrared wavelengths for gold, silver, and aluminum. The dielectric 
constants calculated from that table for the infrared wavelengths appeared to be 
consistent with those obtained in Table V. 

Table VI provides the results of applying the dielectric constants listed in 
Table Y to the single foil case at an energy of 100 MeV. Recall that the total 
intensity should be proportional to the square of the modulus of the difference of the 
dielectric contants of the two media. For the single foil case, the first media was 
assumed to be a vacuum which is characterized by dielectric constant of one, ie e =1 
+ iO. The total intensity ITOT per unit frequency and solid angle in terms of 
charge squared divided by the speed of light is listed for aluminum Al, gold Au, 
and silver Ag. as a function of wavelength (listed in angstroms), and the square of 
the modulus of the difference between the dielectric constant of the metal and of the 
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TABLE V. DIELECTRIC CONSTANT AT OPTICAL WAVELENGTHS 



Metal 


Wavelength 

(angstroms) 


n 


K 


el 


e2 


A1 


5890 


1.44 


3.69 


11.5 


10.6 


Au 


4410 


1.18 


1.56 


1.04 


3.70 


Au 


5890 


0.47 


6.02 


36.0 


5.66 


Ag 


5000 


0.17 


17.1 


292 


5.88 


Ag 


5890 


0.18 


20.6 


424 


7.28 



o 

vacuum (e 9 — e^) . The dielectric constant of the metal for the wavelength listed is 
the same as that listed for the same wavelength in Table V. 



TABLE VI. SINGLE FOIL DIELECTRIC CONSTANT VS. TOTAL 





TRANSITION 


RADIATION INTENSITY AT 100 MEV. 


Metal 


Wavelength 


( 

( £ rV 


ITOT 




(angstroms) 




rtf') 

e /c 
dudD_ 




A1 


5890 


138 


881 


Au 


4410 


8.08 


684 


Au 


5890 


1250 


S45 


A g 


5000 


84700 


964 


Ag 


5890 


179000 


970 



The data in Table VI clearly indicates that the total intensity does indeed 
increase whenever the square of the modulus of the difference in dielectric constants 
is increased. Figures 29—31 depict the intensity distributions of aluminum, gold, 
and silver at a wavelength of 5S90 angstroms and energy of 100 MeV. The data in 
each plot is normalized to the maximum value, represented by the top of the 
z— axis. The data extends five times the predicted angle of maximum intensity in the 
x— y plane. Since the data in each case is plotted on the same relative scale, Figures 
29—31 appear nearly identical. 
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Figure 29. Tlie single foil transition radiation intensity distribution for aluminum 

at 100 MeV at a wavelength of 5890 X. The data extends to ± 5/7 radians in the 
X— V plane. The top of the Z— axis corresponds to the maximum intensity of the 
data. 
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Figure 30. The single foil transition radiation intensity distribution for gold at 100 

MeV at a wavelength of 5890 X. The data extends to ± 5/7 radians in the X— Y 
plane. The top of the Z— axis corresponds to the maximum intensity of the data. 
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Figure 31. The single foil transition radiation intensity distribution for silver at 100 

MeV at a wavelength of 5S90 X. Note that this plot is nearly identical to the plots 
for aluminum and gold (Figures 29, 30). Dielectric permittivity does not affect the 
angular distribution of transition radiation. It does have a significant effect on the 
intensity. The data extends to ±5/7 radians in the X-Y plane. The top of the 
Z— axis corresponds to the maximum intensity of the data. 
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IV. COHERENT TRANSITION RADIATION 



The two— foil interferometer program COHER3 was used to study the effects 
of energy, beam divergence, and dielectric constant on coherent transition radiation 
for comparison with the results obtained from the single foil case. The 
interferometer program permits the study of a finite optical bandwidth of 
wavelength and thus also of frequency. In addition, it promotes a better 

understanding of the concept of coherence length and its importance in determining 
the separation between foils to obtain coherent addition of transition radiation 
from successive foils. 

A. COHERENCE LENGTH 

It is necessary to begin the study of the two— foil interferometer with a 
discussion of coherence length. This is due to the fact that the distance between the 
foils of the interferometer in relation to the coherence length, has a profound effect 
upon the resulting interference pattern. Recall from chapter one that the coherence 
length in a medium is defined 



L 

m 



1 



2 7T | 1— /VfCOS0| 



(4.1) 



Chapter I also described how the forward transition radiation from the back 
of the first foil and the backward transition radiation emitted from the front surface 
of the second foil differ in phase by 
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(4.2) 



$ — P cos 0) — L/L y , 

where L is the distance between the foils and L is the coherence length in a 
vacuum. Obviously, the interference pattern generated by the interferometer is 
highly dependent upon the phase difference and thus the ratio of the distance 
between the foils and the coherence length. For the radiation patterns to add 
coherently, the distance between the foils should be an integer multiple of the 
coherence length in a vacuum for that energy. Wartski [Ref. 5] showed that the 
coherence length on the axis of specular reflection, i.e. 6 = 0, is proportional to the 
square of the energy 

L = 4 Y • (4-3) 

-1 n 

Table VII displays the coherence length on the axis of specular reflection at 
various energies for a wavelength of A = 5890 angstroms. This wavelength was 
chosen because it was the only wavelength in the visible spectrum for which the 
optical constants for determining the dielectric constant were available for gold, 
silver, and aluminum. The default value of the two-foil interferometer program for 
the distance between the foils is 1 .20 centimeters. 

Figures 32 and 33 demonstrate the impact of coherence length on the 
interference pattern from a two— foil interferometer. For both cases the energy is 30 
MeV and the data in the X-Y plane is spread over ± 3/7 radians from the Z-axis. 
The top of the Z-axis corresponds to the maximum value of the data. In Figure 32 
the distance between the foils is 1.20 centimeters while in Figure 33 the distance is 
equal to the coherence length of .034 centimeters. All angular information available 
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Figure 32. Two-foil interferometer interference pattern at 30 Mev for a clear front 
foil and gold back foil. Distance between foils is 1.20 cm. while coherence length is 

.034 cm. at 30 MeV for a wavelength of 5890 X. 
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in Figure 33, where the distance between the foils is equal to the coherence length, is 
lost in Figure 32. Not suprisingly the intensities for each case were completely 
different. A numerical comparison of the intensities would not be meaningful since 
one case is coherent and the other is not. 



TABLE VII. COHERENCE LENGTH VS. ENERGY FOR A = 5890 X 



Energy 


Coherence 


(MeV) 


(cm) 


10 


.004 


20 


.015 


30 


.034 


40 


.059 


50 


.092 


60 


.131 


70 


.178 


80 


.233 


90 


.294 


100 


.363 


150 


.813 


200 


1.44 



B. ENERGY 

The discussion on coherence length demonstrated the effect of energy on the 
considerations for determining foil separation as a function of energy and coherence 
length. Table VIII demonstrates the effect of energy on the angle at which 
maximum intensity occurs, the parallel and perpendicular components of intensity 
and the total intensity for comparison with the effects of the single foil case. The 
intensities listed are for the maximum of each component and therefore should not 
be summed to obtain the maximum total intensity, ITOT. The angle at which the 
maximum intensity occured is much greater than the 1/7 angle predicted for the 
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TABLE VIII. ENERGY EFFECTS ON INTERFEROMETER 
TRANSITION RADIATION 



Energy 


1/7 


Actual 0 ^ 




ITOT 


IPAR 


IPERP 


(MeV) 


(radians) 


(radians) 




2, i 

e /c 




2/ i 

e / c 




2, 

e / c 




duxID^ 




dodfl 




_dodD_ 


20 


.0249 


.0730 


31.92 


31.92 


.0008 


60 


.0084 


.0536 


900.3 


900.3 




1615 


100 


.0051 


.0277 


2046 


2046 


.8868 



single foil case. A comparison with the intensities in Table II of Chapter III 
indicates that the intensities for the interferometer case are much more sensitive to 
increases in energy than for the single foil case. The fact that the intensity of TR 
can be increased through the use of a stack of plates is well known and has proven 
to be a useful tool in studying transition radiation [Ref. 2-4], Ginsburg and 
Tsytovich [Ref. 1] and Garibian [Ref. 4] showed that a condition of resonance occurs 
when the distance between successive foils is such that the transition radiation fields 
from the foils add coherently. This condition is met when the distance between the 
foils is an integer multiple of the coherence length. If the distance between the foils 
does not meet this criteria, the intensity is then a function of the number of 
boundaries per unit length. For this study the separation distance between the foils 
was set equal to the vacuum coherence length for the energy listed. 

Figures 34—36 depict coherent TR patterns at 60 and 100 MeV. The data is 
spread over an angular distribution of five times the predicted angle of maximum 
intensity for the single foil case. The top of the Z— axis corresponds to the 
maximimum value of the total intensity. 
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Figure 34. Coherent TR interference pattern at 60 MeV for a wavelength of 5890 X. 
Separation distance between foils was 1.31 mm. The Z— axis measures relative 
intensity. The top of the Z— axis corresponds to the maximum intensity in units of 
charge squared divided by the speed of light per unit frequency and solid angle. 
Angular information is depicted in the X— Y plane. The data extends to five times 
the predicted angle of maximum intensity (in radians) for the single foil case. 
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Figure 35. Another view of coherent TR interference pattern at 60 MeV for a 

wavelength of 5890 X. Separation distance between foils was 1.31 mm. The Z— axis 
measures relative intensity. The top of the Z-axis corresponds to the maximum 
intensity in units of charge squared divided by the speed of light per unit frequency 
and solid angle. Angular information is depicted in the X-Y plane. The data 
extends to five times the predicted angle of maximum intensity (in radians) for the 
single foil case. 
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Figure 36. Coherent TR interference pattern at 100 MeV for a wavelength of 5890 

X. Separation distance between foils was 3.63 mm. The Z— axis measures relative 
intensity. The top of the Z-axis corresponds to the maximum intensity in units of 
charge squared divided by the speed of light per unit frequency and solid angle. 
Angular information is depicted in the X— Y plane. The data extends to five times 
the predicted angle of maximum intensity (in radians) for the single foil case. 
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The characteristics of the surface plot are quite different than that of the 
single foil case. A series of fringes surround the central cone of maximum intensity. 
The positions of these relative maximum and minimum intensities are described by 
Wartski [Ref. 5] in terms of an interference order p which is defined 

p = -^(1— PcosO) = > (4-4) 



where L is the coherence length in a vacuum and L is the distance between the 
foils. For the relativistic case, this may be approximated as 

P = 5^7 -2 +« 2 ) . (-*'5) 



Local minima. 0 , occur whenever p is an integer, k. Local maximum 
intensities, 6 , occur whenever p = k±l/2. The order of interference at the center, 
that is on the axis of specular reflection is defined to be 



P 



o 




(4.6) 



The angle for which the local maxima and minima intensities occur can be 
described by 



7 M,m 



2A , n 
T^P-P 



o' 



for p = k±l/2 and p = k, respectively. 



(4.7) 
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The intensites for the coherent interferometer case increase with an increase in 
energy. The intensities for the incoherent case depend upon the number of 
boundaries per unit length as well as upon intensity. The maximum intensity 
occurs at the inner most fringe which is located at a much greater angle than that 
predicted for the single foil case. The locations of maxima and minima are 
described in terms of an order of interference p, wavelength A, and distance 
between foils L. The order of interference was shown to be inversely proportional 
to the square of the Lorentz factor 7. 

C. BEAM DIVERGENCE 

The effects of beam divergence were studied for energies of 10, 40, and 100 
MeV for comparison w r ith beam divergence effects at those energies for the single foil 
case. Although an increase in beam divergence decreases the visibility of the 
fringes, the beam divergence required to w r ipe out all angular information was found 
to be much greater than that required for the single foil case. Coherent transition 
radiation is less susceptible to the effects of beam divergence because of the 
amplified intensities due to resonance. 

Table IX contains a comparison of the total, parallel, and perpendicular 
intensities for a few rms beam divergence angles. In all cases, the rms beam 
divergence angle in the X— Z plane and in the Y— Z plane were set equal to each 
other, i.e. In this case the total beam divergence angle is fia, which is the 

value listed in the table in radians. Energy is listed in MeV, the intensities are per 
unit frequency and solid angle in terms of charge squared divided by the speed of 
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light. As in the single foil case, the values listed for each component of intensity are 
the maximum attained for that component and do not necessarily add up to the 
maximum overall intensity, labelled ITOT. 



TABLE IX. BEAM DIVERGENCE VS COHERENT OPTICAL TRANSITI 
RADIATION INTENSITIES AT A = 5890 X 



Energy Beam Div. 
(MeV) (radians) 



ITOT IPAR IPERP 



[ e 7cl 




'e z /c' 




r z / ^ 

e / c 


dodll 

v. J 











Actual 0 

m 

(radians) 



10 


1.4 IE— 2 


23.27 


23.01 


.7195 


.2654 


10 


2.83E-2 


21.87 


20.79 


2.685 


.2654 


10 


4.24E-2 


20.72 


18.14 


5.650 


.2627 


10 


5.66E-2 


20.72 


15.76 


9.356 


.2581 


10 


8.49E-2 


26.75 


12.50 


18.15 


.0000 


40 


7.07E-2 


308.2 


293.6 


37.07 


.0688 


40 


8.49E-2 


300.1 


278.7 


51.94 


.0695 


40 


1.41 E — 1 


2587 


69.93 


2517 


.0000 


100 


2.83E-3 


1877 


1788 


224.9 


.0277 


100 


4.24E-3 


1768 


1560 


472.7 


.0274 


100 


7.07E-3 


1889 


1186 


1145 


.0270 



An increase in beam divergence resulted in a corresponding increase in the 
perpendicular component of intensity and decrease in the parallel component of 
intensity. Also the total intensity tended to decrease until the perpendicular 
component of intensity became comparable to or greater than the parallel 
component. The maximum intensity occurred on the axis of specular reflection 
when the perpendicular component of intensity dominated the intensity 
distribution. The physics behind this was explained for the single foil beam 
divergence in Chapter III. 
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The effect of beam divergence on the coherent transition radiation interference 
pattern is depicted in Figures 37-43. For all the figures, the angular information 
depicted in the X— Y plane extends to five times the predicted angle of maximum 
intensity in radians for the single foil case at the same energy. Intensity is 
measured along the Z— axis. The top of the Z— axis corresponds to the maximum 
intensity in units of charge squared divided by the speed of light per unit frequency 
and solid angle. 

The apparent inconsistency that appears in these plots is due to a forced 
asymmetry built into the program. When a beam of particles encounters a 
boundary at an oblique angle, the resulting transition radiation is assymetric in that 
one side of the pattern will be more intense than the other. This asymetry was 
simulated in the two-dimensional program by changing the sign of the observation 
angle on opposite sides of the axis of specular reflection. That technique was carried 
on into the three-dimensional simulation. Obviously, a refinement of the method is 
needed, but the solution is not a simple one. The inconsistency could be removed 
by neglecting the asymetry introduced by the beam encountering the boundary at 
angle. 

Figures 37-40 depict the effect of beam divergence on the fringe pattern at an 
energy of 10 MeV. As the beam divergence is increased, the fringes become less 
visible. Note that the central cone becomes more shallow as well as the 
perpendicular component of intensity increases. All the angular information is lost 
when the beam divergence angle becomes sufficiently large for the perpendicular 
component of intensity to dominate the intensity distribution. However, the beam 
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Figure 37. Coherent transition radiation at 10 MeV for an rms beam divergence 
angle of 1.41 E — 2 radians. The apparent inconsistency is due to a crude 
approximation of the asymmetry arising from the charged particles encountering the 
boundary at an oblique angle. The angular information depicted in the X-Y plane 
extends to five times the predicted angle of maximum intensity in radians for the 
single foil case at the same energy. Intensity is measured along the Z-axis. The top 
of the Z— axis corresponds to the maximum intensity in units of charge squared 
divided by the speed of light per unit frequency and solid angle. 
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Figure 38. Coherent transition radiation at 10 MeV for an rms beam divergence 
angle of 2.83E-2 radians. Note the relative amplitudes of the outer fringes is much 
less than those of Figure 37. The angular information depicted in the X-Y plane 
extends to five times the predicted angle of maximum intensity in radians for the 
single foil case at the same energy. Intensity is measured along the Z-axis. The top 
of the Z-axis corresponds to the maximum intensity in units of charge squared 
divided by the speed of light per unit frequency and solid angle. 
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Figure 39. Coherent transition radiation at 10 MeV for an rms beam divergence 
angle of 5.66E— 2 radians. Note that the outer fringes are no longer visible. 
Compare with Figures 37 and 38. Note also the general speading and smoothing of 
the intensity distribution. The angular information depicted in the X— Y plane 
extends to five times the predicted angle of maximum intensity in radians for the 
single foil case at the same energy. Intensity is measured along the Z— axis. The top 
of the Z— axis corresponds to the maximum intensity in units of charge squared 
divided by the speed of light per unit frequency and solid angle. 
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Figure 40. Coherent transition radiation at 10 MeV for an rms beam divergence 
angle of 8.49E— 2 radians. The beam divergence has washed out all angular 
information. However a beam divergence angle of little more than 2.50E— 2 radians 
accomplished the same effect for the single foil case. The angular information 
depicted in the X— Y plane extends to five times the predicted angle of maximum 
intensity in radians for the single foil case at the same energy. Intensity is 
measured along the Z— axis. The top of the Z— axis corresponds to the maximum 
intensity in units of charge squared divided by the speed of light per unit frequency 
and solid angle. 
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divergence angle at which this occurs is significantly higher than that for the single 
foil case. Figure 40 depicts the intensity distribution when the beam divergence 
angle is sufficiently large enough for this phenomenon to occur. 

Figures 41 and 42 depict the effect of beam divergence on the interference 
pattern at 40 MeV. Here again, the fringe parrern becomes less visible as the beam 
divergence is increased. Finally, Figure 43 depicts the case where the beam 
divergence is nearly large enough for the angular information to be completely 
washed out. The beam divergence angle of 7.07E-3 radians required for this is well 
over twice that required of the 3.00E— 3 radians for the single foil case at the same 
energy of 100 MeV. 

In general, beam divergence has the same effect on the components of intensity 
an for the coherent transition radiation case as for single foil transition radiation. 
As in the single foil case, the effects of beam divergence are amplified at higher 
energies. An increase in beam divergence decreases the visibility of the fringes. 
However, the effect of resonance causes coherent transition radiation to be less 
susceptible to the loss of angular information due to beam divergence. 

D. DIELECTRIC CONSTANTS 

The intensity of transition radiation emitted when a charged particle transits 
from one medium to another is dependent upon the dielectric constant of each 
media. It was shown in Chapter III that the intensity is directly proportional to the 
square of the modulus of the difference between the dielectric constants of the two 
media. As in the single foil case, a vacuum to medium transition was assumed. The 
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dielcctiic constants for gold Au, silver Ag, and aluminum A1 at a wavelength of 
5890 X were used to study the effect on coherent transition radiation at 100 MeV. 




Figure 41. Coherent transition radiation at 40 MeV for an rms beam divergence 
angle of 7.07E-3 radians. The angular information depicted in the X-Y plane 
extends to five times the predicted angle of maximum intensity in radians for the 
single foil ease at the same energy. Intensity is measured along the Z— axis. The top 
of the Z— axis corresponds to the maximum intensity in units of charge squared 
divided by the speed of light per unit frequency and solid angle. 
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Figure 42. Coherent transition radiation at 40 MeV for an rms beam divergence 
angle of 8.49E— 3 radians. Note the decrease of the outer fringe amplitudes as 
compared to those in Figure 41. The angular information depicted in the X— Y 
plane extends to five times the predicted angle of maximum intensity in radians for 
the single foil case at the same energy. Intensity is measured along the Z— axis. The 
top of the Z— axis corresponds to the maximum intensity in units of charge squared 
divided by the speed of light per unit frequency and solid angle. 
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Figure 43. Coherent, transition radiation at 100 MeV for an rms beam divergence 
angle of 7.07E— 3 radians. The beam divergence angle is nearly great enough to 
wash out all angular information. However, a beam divergence angle of only 
3.00E— 3 radians was able to produce the same effect for the single foil case. The 
resonance effects in coherent transition radiation reduces the susceptibility to loss of 
angular information caused by beam divergence. The angular information in the 
X-Y plane extends to five times the predicted angle of maximum intensity in 
radians for the single foil case at the same energy. Intensity is measured along the 
Z— axis. The top of the Z— axis corresponds to the maximum intensity in units of 
charge squared divided by the speed of light per unit frequency and solid angle. 
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The results of applying the dielectric constants found in Table V of Chapter 
III for the parameters discussed here are presented in Table X for comparison with 
the results of the single foil case. The total intensity ITOT is presented in units of 
charge squared divided by the speed of light per unit frequency and solid angle. 

TABLE X. DIELECTRIC CONSTANT VS. INTENSITY FOR 
COHERENT TRANSITION RADIATION AT 100 MEV AND A=5890 X 



Metal 


( e l 


ITOT 




fe 2 /c 






dudD 


Aluminum A1 


138 


1757 


Gold Au 


1250 


1964 


Silver Ag 


179000 


2256 



The intensities for the coherent transition radiation case are much greater 
than that for the single foil case. This effect is again attributable to resonance. It is 
clear that the intensity for the coherent case is indeed proportional to the square of 
the difference between the dielectric constants between thea two media. Figure 44 
depicts the coherent transition radiation interference pattern for silver at 100 MeY 
over an angular distribution of ± 5/q radians from the Z— axis. The data is 
normalized so that the top of the Z— axis corresponds to the maximum value of 
intensity. 

E. FINITE BANDWIDTH 

Wartski [Ref. 5] found that variations in wavelength, and thus frequency, 
affected the visibility of the interference pattern fringes from the interferometer. 
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The pattern fringes became less visible as the bandwidth of wavelength increased. 
Wartski showed that the relation between intensity per unit solid angle and a 
bandwidth of wavelength could be expressed 



where aA is the bandwith of wavelength, L is the distance between foils, and p is 
the order of interference. An expression for the intensities of the maxima and 
minima may be found by substituting p=k±l /2 or p=I\ respectively. 



bandwidth resulted in very slight fluctuations of intensity. There also appeared to 
be a slight shift in the fringe patterns for the 30 MeV case. 

Figures 45-48 depict the effect of two percent and fifteen percent wavelength 
bandwidths at 30 and 150 MeV. The data is normalized so that the top of the 
Z— axis corresponds to the maximum intensity. The X— Y plane is used to plot the 
angular information in terms of the predicted angle of maximum intensity, 1 / 7 , for 
the single foil case. The data extends to ± 5/7 on the X and Y axes. A close 
comparison of the fringe patterns at 30 MeV for the four percent bandwidth in 




(4.8) 



For this study, the dielectric constants of gold at wavelengths of 4410 X and 
5890 X was used. The coherent transition radiation distribution was studied at 30 
MeV for a single wavelength of 4410 X, and for a two percent (88 X) and fifteen 
percent (66 X) bandwidth centered on this wavelength. A similar study was done at 
150 MeV around a wavelength of 5890 X. It was found that an increase in 
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Figure 44. Coherent TR interference pattern for silver at 100 MeV and a 

wavelength of 5890 X. The patterns for silver and gold at the same eneigy and 
wavelength are nearly identical despite the difference in delectric permittivity. 

The reason for this is that the data in each case is plotted on a relative scale. The 
intensity is normalized so that the top of the Z-axis corresponnds to the maximim 
intensity. The data in the X— Y plane extends to five times the predicted angle of 
maximum intensity for the single foil case at 100 MeV. 
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Figure 45. Coherent TR interference pattern at 30 MeV for a four percent 

wavelength bandwidth centered at 4410 X. The data is normalized so that the top 
of the Z— axis corresponds to the maximum intensity. The X— Y plane is used to 
plot the angular information in terms of the predicted angle of maximum intensity, 
1 / 7 , for the single foil case. The data extends to ±5 / 7 on the X and Y axes. 
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Figure 46. Coherent Til interference pattern at 30 MeV for a fifteen percent 

wavelength bandwidth centered at 4410 X. A close comparison of the fringe 
patterns here with those in Figure 45 reveals that the fringe patterns become 
slightly less visible as the bandwidth is increased. The data is normalized so that 
the top of the Z-axis corresponds to the maximum intensity. The X-Y plane is 
used to plot the angular information in terms of the predicted angle of maximum 
intensity, l/% for the single foil case. The data extends to ± 5/7 on the X and Y 
axes. 
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Figure 47. Coherent TR interference pattern at 150 MeV for a four percent 

wavelength bandwidth centered at 5890 X. The data is normalized so that the top 
of the Z-axis corresponds to the maximum intensity. The X— Y plane is used to 
plot the angular information in terms of the predicted angle of maximum intensity, 
1/7, for the single foil case. The data extends to ± 5/7 on the X and Y axes. 
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Figure 48. Coherent TR interference pattern at 150 MeV for a fifteen percent 

wavelength bandwidth centered at 5890 X. This pattern is indistinguishible from 
that of the four percent bandwidth presented in Figure 47. The data is normalized 
so that the top of the Z-axis corresponds to the maximum intensity. The X— Y lane 
is used to plot the angular information in terms of the predicted angle of maximum 
intensity, 1 / 7 , for the single foil case. The data extends to ± 5/7 on the X and Y 
axes. 
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Figure 45 and the fifteen percent bandwidth in Figure 46 reveals that the outer 
fringes become slightly less visible as the bandwidth is increased. However, no such 
difference could be found for the same case at 150 MeV (Figures 47,48). 

In summary, a finite bandwidth of wavelength, and thus frequenccy, results in 
slight variations of intensity. An increase in bandwidth produces a small increase in 
intensity. Angular information such as fringe pattern is relatively unaffected except 
for large bandwidths. 
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RECOMMENDATIONS AND CONCLUSIONS 



The programs used in this study provide an effective tool for studying the 
characteristics of optical transition radiation to aid in interpretation of experiments. 
They provide an inexpensive method of obtaining a thorough understanding of 
optical transition radiation and its potential use for charged particle beam 
diagnostics. Such an understanding has already been promoted through the results 
of this study. 

The nature of the energy effects on angular distribution and polarization was 
thoroughly explored. Other factors affecting the intensity and distribution of 
transition radiation were also investigated It was shown that beam divergence 
affects the angular distribution of intensity, particularly at higher energies. 
However, the effects of beam divergence can be more sensitively investigated 
through the use of a stack of foils spaced at intervals comparable with the coherence 
length. Metals characterized by dielectric prermittivities much greater than one can 
be used to enhance the intensity of transition radiation. The finite bandwidths of 
the observed wavelength have little effect on the distribution of transition radiation 
unless the bandwidth is excessively large. 

Further improvements can be made on the programs used in this study. For 
the single foil and two— foil interferometer TR programs, the asymetry representing 
the effect of a charged particle encountering a boundary at an oblique angle requires 
refinement. A cosine dependence was discussed in which the asymmetry could 
beintroduced at one point and then "rotated" around the axis. The asymmetry 
factor could then be multiplied by the cosine of the angle of rotation. Some means 
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of labelling or indicating units of measure on the axes for the plots in the SON OF 
SURF3D program is also recommended. The SON OF SURF3D program should 
also be fully developed for commercial use. 

Finally, software should be developed to measure and present actual data in 
three dimensions for comparison with the theoretical data presented here. A 
program compatible with SON OF SURF3D under development for that purpose 
had to be abandoned due to the time constraints on this study. Such a program 
would be a necessary step in actually applying the optical transition radiation as a 
diagnostic for charged particle beams. 
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APPENDIX A 



USER'S MANUAL 

This appendix is intended to aid the operator in using the transition radiation 
programs TR3 and COHER3 and the three— dimensional graphics program SON OF 
SURF3D. While the content was written so as to be readily understandable to the 
inexperienced computer user, a basic understanding of MS-DOS procedures is 
assumed. This appendix is written in two parts. The two transition radiation 
programs, being very similar in use are treated as one unit. A description for the 
use of each was taken up to the point of invoking the use of the graphics program. 
The final section is written so that the user would be able to begin with a data file 
from either transition radiation program, plot it out in three dimensions, and vary 
all the plotting parameters in any combination so as to obtain the most satisfactory 
plot for the case at hand. 

A. THE TRANSITION RADIATION PROGRAMS 

The transition radiation programs were written in Fortran and compiled using 
the Ryan— McFarland fortran compiler. Both programs should be installed on the 
hard disk under the directory C:\FORT>. If not previously compiled, the 
programs may be compiled from within the C:\FORT> directory using the 
command RMFORT TR3 /N or RMFORT COHER3 /N. This command causes 
and object file to be formed which is used by the PLINK command, PLINK86 FI 
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FR3 or P LINK 86 FI COHER3 at the C:\FO RT> prompt, to form an execute file. 
The programs are then ready for use by typing either TR3 or COHER3 for the 
program desired. 

1. Overview 

Program execution begins upon typing either TR3 or COHER3 at the 
C:\FORT> prompt. The user is then faced with the following messages in 
succession: 

ENTER NAME OF 3-D GRAPHICS FILE 

EXAMPLE: GRAPH.DAT 

[entrv] 

ENTER NAME OF CORRESPONDING PARAMETER FILE 

EXAMPLE: PRMTRS.DAT 

[entry] 

WHAT IS THE DESIRED BEAM ENERGY IN MEV? 

PLEASE INCLUDE A DECIMAL PLACE. 

[entry] 

Once the above entries are made, a menu list of current parameters 
appears. This menu list for the COHER3 program differs a bit from that of the 
TR3 program as should be expected. The menu list for the TR3 and the COHER3 
programs are depicted in Figures 49 and 50, respectively. The user has the option 
of accepting the current default parameters or changing one, any combination, or all 
the parameters. In order to change any parameter, the user need only type the 
appropriate number, press return/enter, and follow the directions. In all cases, after 
the new parameter is entered, the parameter change menu reappears indicating the 
new parameters. Once all the parameters are to the user's satisfaction, the user 
accepts the current parameters by entering the appropriate number as indicated in 
order to continue with the program. 
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The last entry required of the user provides the user with a choice of 
polarization. The parallel component of intensity only, the perpendicular 
component of intensity only, or the sum of the two for the total, non— polarized 
intensity may be plotted. The decided choice is entered by typing the appropriate 
number and pressing return/enter. 

CHOOSE DESIRED POLARIZATION. 

ENTER APPROPRIATE INTEGER FOR DESIRED POLARIZATION 

1 NO POLARIZATION— TOTAL INTENSITY 

2 PARALLEL. 

3 PERPENDICULAR 

Once the choice of polarization is entered, the program begins 
calculations that take approximately seven to ten minutes (assuming a math 
co— processor is installed). Progress is indicated by the program through displaying 
the current record number over which the calculations are being done. There are a 
total of 64 records corresponding to rows of data. The data is organized into 64 
rows (records) and 64 columns (fields) for use with the graphics program. All data 
for the graphics file is normalized after the intensities for all the data points are 
calculated. The output includes the graphics file of normalized data., a 
corresponding parameter file, a file containing the non— normalized intensity data in 
terms of charge squared divided by the speed of light (per unit frequency and solid 
angle), and a scratch file used for interim calculations. 

2. Required entries 

The program is executed by typing TR3 or COHER3 at the C:\FORT> 
prompt. The following message then appears upon the screen: 

ENTER NAME OF 3-D GRAPHICS FILE 

EXAMPLE: GRAPH.DAT 
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The required entry here is the name of the actual file that will be used 
by the SON OF SURF3D graphics program. It is suggested that some sort of 
mnemonic device be used for key characteristics of the file for ease in later 
identification. Table XI demonstrates the standard format used in conjunction with 
this thesis. 



TABLE XI. GRAPHICS FILE LABELLING SCHEME 

AXXXENXX.DAT A Either T for TR3 single— foil program or 

C for COHER3 interferometer program. 

.DAT Extension required by graphics program. 

EXXX Indicated energy in MeV. 

XXX Specified key characteristics. 

3 Standard defaults, ie T3E10.DAT 

AN RMS Beam angle altered, N indicated 
variation number for that energy, ie 
CA2E30.DAT for second variant at 30 MeV. 
F scattering foil present, ie TFE50.DAT 

DXX Dielectric constant altered, ie CDA1 for 
dielectric constant of Aluminum. 

NG Data calculated over angular distance of 
N times 1/7 from axis, ie T5GE70.DAT 
N Data normalized to maximum value for 

total intensity at 100 MeV. ie CNE90.DAT 
PERP Perpendicular component of intensity 

only selected for plots, ie TPERPE5.DAT 
PAR Parallel component of intensity only 
selected for plots, ie CPARE15.DAT 
B Finite bandwidth, wavelengths altered 

COHER3 program only, ie CBE80.DAT 
TH Foil thickness altered, COIIER3 program 
only. ieCTHE50.DAT 

DL Distance between foils altered, COHER3 
program only, ie CDLE90.DAT. 



Combinations of key characteristics are also possible. For example, 
TF0GNE8O.DAT indicates the presence of a scattering foil, calcualtions over an 
angular distance of ±(5/7) from the axis, and data normalized to 100 MeV. 
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Upon entering the name of the graphics file, the operator is requested to 
enter the name of the associated parameters file. 

ENTER NAME OF CORRESPONDING PARAMETER FILE 
EXAMPLE: GRAPHP.DAT 

The name of the parameter file should be similar enough to the name of 
the graphics file to avoid confusion at a later date. For this thesis a P for 
parameters was tagged on to the name of the graphics file. For example, 
C5GE60P.DAT or C5GP60.DAT for the parameter file corresponding to the 
C5GE60.DAT graphics file. 

One word of caution. The computer doesn't read beyond the first nine 
letters of a file name. The file name plus three— letter extension should never exceed 
twelve letters. The programs will run using filenames entered in excess of nine 
characters. However, the files in which the only difference in name occurs beyond 
the ninth character will not appear in the directories and the data will be lost. It 
may occassionally become necessary to abbreviate the name of the graphics and 
parameter file to prevent this from happening. 

Once the name of the parameter file has been entered, the user is 
requested to enter the energy in MeV. Since the analytic solutions were obtained 
through the use of small angle approximations for angles represented by the inverse 
of the Lorentz factor, analysis at lower energies must be considered less exact. 
Energies at which the inverse of the Lorentz factor is greater than 0.1, that 
isenergies less than about 5 MeV should be treated as rough approximations. 
Angular information, which is important for purposes of beam diagnostics, is lost 
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above about 200 MeV without a polarizing filter. No such filter is built into these 
programs therefore their analysis above 200 MeV is of little value for beam 
diagnostics. The energy entry should include a decimal point. 

Once the output files have been named and the desired energy entered, a 
menu for changing data parameters appears on the screen. Figures 49 and 50 depict 
the menus for the TR3 single— foil and the COHER3 interferometer programs 
respectively. 



The current values are: 

1 RMS beam angle is 0.001 radians/sqrt(2). 

2 IFOIL is 0. 

A value of one indicates that a scattering foil 
is present. A value of zero indicates that there 
is a target foil only. 

3 The dielectric constant of the medium is 

0.618 5.4 1 i. 

This program assumes a vacuum to medium transition. 

4 Angular measurement is over 0.05 radians. 

5 Normalization is to the maximum total 
intensity for the energy entered for this run. 

The option is to normalize to the maximum total 
intensity for an energy of 100 MeV. 

In all cases, parallel intensity will be 
normalized to its own maximum value for the 
energy entered for this run. 

6 Accept current parameters. 

Enter the number without the decimal for the value 
you wish to change. Enter 6 for no changes. 

Figure 49. Parameter change menu for the single foil TR3 program 
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THE CURRENT PARAMETERS ARE: 



1 RMS beam angle in the X— Z plane is .000210 radians 
RMS beam angle in the Y— Z plane is .000210 radians 

2 IFOIL is 0. A value of 1 indicates 

that a scattering foil is present. A value of 
zero indicate absence of scattering foil. 

3 The optical constants for the back metalic foil 
are .618 + 5.47i. 

4 The optical constants for the front clear foil 
are 1.48 + O.OOi. 

5 Angular measurement is over 0.05 radians 

6 Normalization is to the maximum total 
intensity for the enery entered for this run. 

The option is to normalize to the maximum total 
intensity for an energy of 100 MeV. 

7 Minimum wavelength is 4500 angstroms. 

Maximum wavelength is 4500 angstroms. 

8 Foil thickness is .005 centimeters. 

9 Distance between foils is 1.20 centemeters. 

10 Accept current parameters. 

Enter the number without the decimal for the value 
you wish to change. Enter 10 to accept parameters. 

Figure 50. Parameter change menu for interferometer COHER3 program 



The user has the option of accepting the current default parameters or 
changing one, any combination, or all the parameters. In order to change any 
parameter, the user must first type the appropriate number, then press 
return/enter. The user will then be directed to enter the new value. Each of these 
parameters will now be discussed in detail. 

Any parameter may be changed simply by entering the integer value 
that appears next to the parameter the user wishes to change. The first parameter 
option in both programs is RMS beam angle. Upon entering a 1 to change this 
parameter, the message 'PLEASE ENTER NEW VALUE FOR RMS BEAM 
ANGLE' appears. The appropriate entry is entered in terms of radians. Not all the 
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particles within the beam travel parallel to each other. The particles tend to 
diverge at various angles within the beam. The rms beam angle is the mean of a 
Gaussian distribution of these angles of divergence. For the TR3 program it is 
assumed that the distribution of divergence in the X-Z plane (a ,) is the same as in 
the Y— Z plane (<r ) such that a =a . In that case the total RMS beam angle of 

*/ J 

divergence is equal to -J2a. The value entered for the TR3 program should be in 
terms of radians/V?. The COHER3 program does not make that assumption and 
allows separate values to be entered for a and o The values entered for the 

x y 

COHER3 program should be in radians. 

Option 2 in both programs allows for the inclusion of a separate 
scattering foil. The message 'ENTER NEW VALUE FOR IFOIL' appears when 
this option is selected. The only valid entries for this option are 0 and 1. As 
indicated, a value of 1 here indicates the presence of the scattering foil. The 
presence of the scattering foil slightly affects the distribution of the parallel 
component of TR intensity. It does not affect the perpendicular component of 
intensity. A value of 0 here indicates that a scattering foil is not present. 

The ability to change the dielectric constants of the mediums is offered 
in option 3 for the single foil TR3 program and by options 3 and 4 of the 
interferometer COHER3 program. Upon selecting any of these options, the user is 
first asked to enter the value of the real component, then asked to enter the 
imaginary component. In effect, this allows the user to change the material of the 
target foil. Option 4 of the COHER3 program allows the user to substitute another 
material for the clear front foil figured into the program. Values for dielectric 
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constants may be found in the CRC Handbook of Chemistry and Physics [Ref 9]. 
Both the real and imaginary parts of the dielectric constant should be entered as 
real numbers. 

Option 4 of the TR3 program and 5 of the COHER3 program allow the 
user to change the anglular range over which the data is calculated. Two different 
means are allowed for this option. The predicted angle of maximum intensty is that 
of the inverse of the Lorentz factor in radians. The user has the option of entering a 
multiple of this angle of entering a value for this angle in radians directly. Upon 
selection of this option, the following menu appears. 

Choose option 

1 Angle that is N*(l /Lorentz factor) 

2 Enter angle in radians manually. 

Choosing option 1 results in the message 'Enter N, for N*(l/Lorentz 
factor)'. The preferred entry is an integer. Choosing this option causes calculations 
of intensity to be over an integer multiple of the predicted angle of maximum 
intensity. This option aids the comparison of energy distributions over a wide range 
of energies. Choosing option 2 results in the message 'Enter desired angle in 
radians'. Obviously, the desired number is a real value of radians for the desired 
angle. Using this option to choose the same angle for a range of energies aids in 
examining the angular dependence of the intensity distribution on energy. 

The ability to normalize data to the maximum intensity value attained 
at 100 MeV is offered by options 5 and 6 of the TR3 and COIIER3 programs 
respectively. This option allows for ease of study of relative intensities as a function 
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of energy. Choosing this option is irreversible. Once chosen, there is no way of 
going back to normalizing the data to its own maximum intensity for that energy. 
Upon entering the appropriate number to choose this option, the parameter change 
menu appears with the message 'Normalization will be to 100 MeV.' in place of this 
option. 

Option 7 of the COHER3 program allows the user to change the range 
of wavelengths over which calculations are made. Selection of this option results in 
the following messages: 

ENTER NEW MINIMUM WAVELENGTH. 

[entrv] 

ENTER NEW MAXIMUM WAVELENGTH. 

[entry] 

The user should enter the minimum and maximum wavelengths in 
angstroms. This option provides for the rare opportunity to study the effects of 
finite bandwidth on the distribution of TR intensity. The traditional development 
of transition radiation theory is built upon an evaluation of integrals that assumes a 
Dirac function for frequency bandwidth. However, a finite frequency bandwidth has 
always been encountered while obtaining real data which may account for the small 
differences between theoretical and actual data. 

Choosing option 8 of the COHER3 program results in the message: 
ENTER NEW VALUE FOR FOIL THICKNESS IN CENTIMETERS. This is for 
the thickness of the front clear foil through which the beam must travel before 
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encountering the back metalic foil. This option provides the opportunity to study 
the effects of front foil thickness on the interference pattern between the forward 
radiation from the front foil and the backward radiation from the rear metallic foil. 

Option 9 of the COHER3 program allows the user to specify the 
distance between foils in centimeters. This option not only provides the user the 
ability to study coherence length, but is anecessity since the interference pattern is 
highly dependent on the distance between foils with respect to the coherence length. 
The convention has been to define coherence length as the distance a particle 
travelsin the time it takes the particle to "see" the phase of its associated 
electromagnetic wave change by one radian [Ref. 5]. This phase difference was 
chosen so that the particle's field and the radiation fields add coherently, hence the 
name "coherence length". 

The programs allow the user to alter any or all data parameters in any 
order or combination. Every time a user chooses an option other than ACCEPT 
CURRENT PARAMETERS a message appears which attempts to clearly specify 
the appropriate entry. Once that entry is made, the parameter change menu 
reflecting the new values reappears on the screen. This allows the user to review all 
the parameter values before continuing on with the program. Once satisfied with 
the current parameters, the user must select ACCEPT CURRENT PARAMETERS 
(option 6 for the TR3 program, option 10 for the COHER3 program) in order for 
the program to continue. 

Having accepted the current parameters, the user is requested to make 
one last selection. The transition radiation intensities are calculated by computing 
the parallel and perpendicular components separately. Total intensity is then the 
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sum of the two components. The user is provided the option of plotting either the 
total intensity or either component of intensity. This allows the user to study the 
effects of each parameter on the components of intensity separately. The user 
selects the desired coomponent (polarization) by entering the appropriate integer as 
specified on the following menu: 

ENTER APPROPRIATE INTEGER FOR DESIRED POLARIZATION 

1 NO POLARIZATION — TOTAL INTENSITY 

2 PARALLEL 

3 PERPENDICULAR 

Once the selection is entered, the program will run its course, indicating 
progress by listing the current record number over which calculations are being 
made, recall that the data is organized into sixtvfour records each containg sixty 
four fields of data for use with the graphics psrogram. The intensities are 
normalized for the graphing purposes after all the other calculations have been 
completed. 

B. THE THREE-DIMENSIONAL GRAPHICS PROGRAM 

The three-dimensional graphics program used in conjunction with the 
transition radiation programs is SON OF SURF3D which was written by Don 
Gilbert of DogStar Software. SON OF SURF3D is an improved version of 
SURF3— D / SURFS7 that is not yet fully developed for public release at the time of 
this writing. The SURF3— D / SURF87 and the improved SON OF SURF3D are 
three dimensional plotting programs that are written in Turbo Pascal and use 
TurboHALO graphics routines [Ref. 10]. Turbo Pascal is a registered trademark of 
Borland International and TurboHALO is a trademark of Media Cybernetics, Inc. 
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and IMSI. The key improvement of SON OF SURF3D over SURF3— D that makes 
it useful for this application is that it is able to plot a three dimensional surface 
from a data file provided by the user. The data file used for this purpose must 
contain the following information in order 

Integer value for the number of y— rows 
Integer value for the number of x— columns 
Real value for the z— axis minimum 
Real value for the z— axis maximum 

Real value for the z— axis center (the intercept with the xy— plane) 

Real z values for each x,y point 

Unfortunately, Don Gilbert of Dogstar Software has indicated that 
development of SON OF SURF3D for public release has stopped. However a copy 
may be obtained for a fee by writing to him at the following address: 

DogStar Software 
PO Box 302 
Boomington, IN 47402 
Telephone (812) 333—5616 

A graphics device and a printer are required for the graphics program. A 
math coprocessor is preferred but not absolutely necessary. SON OF SURF3D 
supports several graphics devices and printers with specific device driver files. The 
appropriate graphics driver device file must be stored in the directory 
C:\T\HALO\ as HALODEF.DEV. The printer device driver file must be renamed 
HALODEF.PRN and placed under the directory C:\T\HALO\ . The files 
SSURF.COM and HALORTP.EXE must be installed under the directory 
C:\T\HALO. All the required files are found in the HALO. ARC file on the disk 
supplied by Dogstar Software. 
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1. USE 



Once installed, the program may be run be typing SSURF from within 
the C:\T\HALO directory. When ready, the program will display the menu 
depicted in Figure 51. The user may choose to plot any of the supplied pre-defined 
functions by pressing the appropriate letter such as C for the Cos Sqr function. 
Choosing any of the pre— defined functions in this way automatically places the 
cursor in the parameters section of the menu. Any or all the parameters may be 
then be altered prior to plotting. The user may choose to plot the program by 
pressing the function key F10 or he may choose another function by pressing Esc 



C' 


) Cos Sqr 


3— D Functions: Z=f(X,Y) 

z=cos( sqr(x)+sqr(y)) 


Ej 


1 Error 


z=erfi ( x, y) 


L) 


Ln Cos 


z=ln(cos(x))— ln(cos(y)) 


Y 


) Your Formula 


z=exp(— sqr(x)/2— sqr(y)/2) 



Plot radius 


= 0.90 


Hatch Increment 


= 0.03 


Delta Increment 


= 0.05 


Hide Back Lines 


= No 


Hide Axes 


= No 


Do X— hatching 


= Yes 


Do Y-hatching 


= No 


Points only 


= No 


Virtual Screen 


= No 



Parameters 

Z— axis Maximum=2.400 
X— axis Maximum=3.000 
V— axis Maximum=3.000 
Z— axis zero =0.000 
X— axis zero =0.000 
Y— axis zero =0.000 
A rotation =3.000 
B rotation =—50.0 
C rotation =—70.0 



' Esc=Exit F2=Print F6=File F10=Draw Arrows=Rotate 



Figure 51. SSURF Main Menu 
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which will place the cursor back into the top part of the menu. Option Y) in the 
top part of the menu is the option of interest in that this option allows the user to 
specify a function or a complete pathname of a data file containing data as described 
above. Choosing this option does not place the cursor in the parameter section of 
the menu. Parameters may not be altered without first obtaining a plot upon 
selection of option Y). Therefore it is recommended that the user first select 
another option as described above, alter the parameters as desired, and press Esc to 
get back into the top portion of the menu before selecting option Y). In order to 
plot the three-dimensional surface of a data file, select option Y) and press 
return/enter. Then type in the complete pathname of the data file and 
pressreturn/enter again. After a few seconds, the word "file" will appear behind the 
pathname indicating that the program is ready to plot the data. The data may then 
be plotted by pressing the function key F10. Flashing indices on the z— axis indicate 
that the plot is finished. The plot may be rotated at will by pressing the arrow 
keys. The user may print the plot by pressing the F2 function key. The Esc key 
will return the display back to the main menu and place the cursor in the 
parameters section. 

The shape of the plot is largely dependent on the plotting parameters 
used to define it. Figure 52 depicts the parameters used for the plots displayed in 
this thesis. 
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Plot radius 


= 0.70 


Hatch Increment 


= 0.05 


Delta Increment 


= 0.05 


Hide Back Lines 


= Yes 


Hide Axes 


= Yes 


Do X— hatching 


= Yes 


Do Y— hatching 


= Yes 


Points only 


= No 


Virtual Screen 


= No 



Parameters 

Z— axis Maximum= 1.000 
X— axis Maximum=40.00 
Y— axis Maximum=40.00 
Z— axis zero =—.500 
X— axis zero =0.000 
Y— axis zero =0.000 
A rotation =3.000 
B rotation =—70.00 
C rotation =-€0.00 



Figure 52. Sample parameter values for plots. 

"Plot radius" refers to the size of the plot in the XY plane relative to the 
extent of the X and Y axes. Valid entries for "Plot radius: range from 0.01 to 1.00. 
"Hatch increment" refers to the spacing between the plotting lines. "Delta 
increment" refers to the increment used for interpolation between data points. 
Valid entries for both "Hatch increment" and "Delta increment" are 0.00 to 1.00. 

The "Hide Back Lines" and the "Hide Axes" options may be used to 
keep the plot relatively clean. The user has the option of selecting either x hatching 
or y hatching or both. Choosing neither produces the same result as "Points only", 
which overides the "X hatching" and "Y hatching" options. The "Virtual Screen" 
option supports the high— resolution capabilities of some printers. A device file 
named HALORVRI.DEV is required for this optioni. Valid entries for all these 
options are either "Yes" or "No". 

"Z— axis maximum" defines the maximum value of the vertical scale, 
regardless of the maximum value of the function or data. "X— axis maximum" and 
"Y— axis maximum" define the maximum horizontal limits of the plot, regardless of 
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the extent of the data. SON OF SURF is capable of plotting as much as a 100 by 
100 array of data centered at the intersection of the three axes. The valid entries 
for these options range from 0.01 to 100.00. 

"Z— axis zero", "X— axis zero", and "Y— axis zero" define the origin of 
respective axis relative to the origin. Selecting something other than zero for any of 
these options physically moves the entire plot. The range of valid entries were not 
listed in the scanty documentation. 

All plots may be rotated about the three axes at will, either by 
specifying the rotation in the parameters section of the menu or by using the arrow 
keys once a plot has been generated. "A rotation" appears to refer to rotation about 
the Z— axis, "B rotation" about the X— axis, and "C rotation" about the Y— axis. 
Although the range of valid entries for these rotations are —360.00 to +360.00, all 
possible rotations can be specified within the range of —180.00 to +180.00. 
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APPENDIX B 



PROGRAM LISTINGS 



A. SINGLE FOIL TR3 PROGRAM 

PROGRAM TR3 

MODIFICATION OF PROGRAM TRADsl, WHICH WAS WRITTEN BY 
D. W. RULE NSWC, CODE MODIFIED SEPT. 1986 TO INCLUDE I_PERP 
modified sept. 1987 to make w(i)= itot=ipar+iperp 

instead of the ratio of iperp to ipar 
MODIFIED BY WIL LONGSTAFF, NPGS, AUG. 1988 
FOR THREE-DIMENSIONAL ANALYSIS 

THIS PROGRAM ASSUMES A VACUUM TO MEDIUM TRANSITION 

CALC. INTENSITY DIST. AS FUNC . OF ANGLE OF TRANSITION RAD. 
ASSUMING GAUSSIAN BEAM SCATTERING ANGLE 

G = LORENTZ FACTOR, S = RMS BEAM ANGLE IN RADIANS /SQRT(2) 

T = ANGLE OF OBSERVATION IN RADIANS 

INTEGER AC, BC, NNW, ROW, COL, XO, YO, COUNT, WC , LW, NC, NW, NI , NWW 
INTEGER ENTRY, CHOICE, POLAR 

REAL THETA (64) , RAD2, RAD, SCALER, ZMIN, ZMAX, ZCEN, TMAX, AMAX, ANG 
REAL W (600) , X(600) , Y (600) , Z (600) , TAU (600) , INT (600) ,MAXY,MAXP, 
+ c(4, 600) ,dx (600) , dy (600) , cy (64 0) , NMBR, PMAX, RMAX 

CHARACTER ANS* 1 , ch * 1 , dat f il * 20, SFILE* 15, DFILE* 15 
COMPLEX AI , EPS, RPAR, RPERP 

DATA PI/3. 141592654/, PSI/0. 7853982/, S/0.001/, IFOIL/0/ 

DATA XN/0 . 618/, XK/5. 47/, IPRINT/0/, INORM/1/, TWOFOIL/1 . 0/ 

DATA XO/ 32 /, YO/ 32/, SCALER/ 620 . 258423/, ZMIN/0 . /, ZMAX/1 . 0/ 

DATA ZCEN/0 . 0/ , ROW/ 64 / , COL/ 64 / , ANG/0 . 05/ , RMAX/0 . 0 / , XDUM/0 . / 
DATA YMAX/0.0/, WMAX/0 . 0/ , PMAX/0 . 0/ , MAXY/0 . 0/ , MAXP/0 . 0 / 

AI = (0 .0, 1 . 0) 

WRITE (6, *) 'ENTER NAME OF 3-D GRAPHICS FILE.' 

WRITE (6, ‘ )' EXAMPLE: GRAPH . DAT ' 

WRITE (6, *) ' ’ 

READ (6, 17) SFILE 

WRITE (6,*) 'ENTER NAME OF CORRESPONDING PARAMETER FILE.’ 

WRITE (6, *) ’EXAMPLE: PRMTRS.DAT’ 

WRITE (6, *) ’ ' 

READ (6, 17) DFILE 
17 FORMAT (A15) 

OPEN ( 3 , F I LE= ' RADOUT . DAT ’ , STATUS= ' UNKNOWN ' ) 

OPEN ( 4 , FILE=DFILE, STATUS= ' UNKNOWN ’ ) 

OPEN (5, FILE=SFILE, FORM= ' FORMATTED ' ) 

OPEN (8 , FILE= ' SCRATCH ’ , STATUS= ' UNKNOWN ’ , ACCESS= 'DIRECT ' , 

+ RECL=256) 

REWIND 2 
REWIND 3 

WRITE (5, 123) ROW, COL, ZMIN, ZMAX, ZCEN 
123 FORMAT (213, 3F8. 4) 
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1 WRITE (6, *) 'WHAT IS THE DESIRED BEAM ENERGY IN MEV? ' 

WRITE (6, *) 'PLEASE INCLUDE A DECIMAL PLACE.' 

WRITE (6, * ) ' VALID ENTRIES ARE 5 MEV - 1000. MEV. ' 

READ (6,*) ENERGY 

IF '{ENERGY . LT . 5 . OR . ENERGY . GT . 1 000 . ) THEN 
GOTO 1 
END IF 

G = 1 .0 + ENERGY/0.511 
TMAX=1 . /G 

2 WRITE ( 6, * ) ' The current values are:' 

WRITE ( 6, * ) ' 1 RMS beam angle is ’,s,' radians/sqrt (2) . ' 

WRITE ( 6, * ) '2 IFOIL is \ifoil,'. ' 

WRITE (6, * ) ' A value of one indicates that a scattering foil' 
WRITE (6, * ) ' is present. A value of zero indicates that there' 
WRITE (6, * ) ' is a target foil only.' 

WRITE ( 6, * ) '3 The dielectric constant of the medium is' 

WRITE (6, *)' ' , XN, ' + ' , XK , ' i . ' 

WRITE (6, *)' This program assumes a vacuum to medium transition.' 
WRITE ( 6, * ) '4 Angular measurement is over’,ANG,' radians.' 

IF (ABS (RMAX-0 . ) . LE . IE- 6) THEN 

WRITE(6,*)'5 Normalization is to the maximum total' 

WRITE ( 6, *)' intensity for the energy entered for this run.' 
WRITE ( 6, * ) ' The option is to normalize to the maximum total' 
WRITE ( 6, * ) ' i ntensi ty for an energy of 100 MeV. ' 

WRITE (6, * ) ' In all cases, parallel intensity will be' 

WRITE (6, *)' normalized to its own maximum value for the' 

WRITE ( 6, *)' energy entered for this run.’ 

ELSE 

WRITE (6, *) 'Normalization will be to 100 MeV.' 

ENDIF 

WRITE (6, * ) ' ’ 

WRITE(6,*)'6 Accept current parameters.' 

WRITE (6, * ) ' ' 

WRITE (6, *) 'Enter the number without the decimal for the value' 
WRITE (6, *)' you wish to change. Enter 6 for no changes.' 

WRITE ( 6, *) ' ' 

READ (6,*) ENTRY 
IF (ABS (ENTRY-6) . LE. IE-6) THEN 
GOTO 3 
ENDIF 

IF (ABS (ENTRY-1) .LE. IE-6) THEN 

WRITE(6,*) 'PLEASE ENTER NEW VALUE FOR RMS BEAM ANGLE.' 

WRITE ( 6, *) ' ’ 

READ (6,*) S 
GOTO 2 
ENDIF 

IF (ABS (ENTRY-2) .LE. IE-6) THEN 
WRITE (6, *) 'ENTER NEW VALUE FOR IFOIL.’ 

WRITE (6, *) ' ' 

READ (6,*) IFOIL 
GOTO 2 
ENDIF 

IF (ABS (ENTRY-3) .LE. IE-6) THEN 

WRITE (6,*) 'Enter the real component: ' 

WRITE (6, * ) ’ ' 

READ (6,*) XN 

WRITE ( 6, *)' Enter the imaginary component without the "i".' 
WRITE(6, *) ' ’ 

READ (6,*) XK 
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GOTO 2 
END IF 

IF (ABS (ENTRY-4) . LE. IE-6) THEN 
WRITE (6, *) 'Choose option.' 

WRITE(6,*)'l Angle that is N*(l/Lorentz factor)’ 

WRITE(6, *)'2 Enter angle in radians manually.' 

WRITE (6, *) ' ' 

READ (6,*) OPTION 

IF (ABS (OPTION-1) .LE . IE-6) THEN 

WRITE ( 6, *)' Enter N, for N*(l/Lorentz factor)' 

WRITE (6, * ) ’ ’ 

READ (6,*) NMBR 
ANG=NMBR/G 

SCALER=31 . / (SIN (ANG) ) 

C WRITE (6, *) 'SCALER IS ', SCALER 

END IF 

IF (ABS (OPTION-2) . LE. IE-6) THEN 

WRITE ( 6, *)’ Enter desired angle in radians.' 

WRITE ( 6, *) ' ' 

READ (6,*) ANG 
SCALER=31 . / (SIN (ANG) ) 

C WRITE (6, *)’ SCALER IS ', SCALER 

ENDIF 
GOTO 2 
ENDIF 

IF (ABS (ENTRY-5) .LE. IE-6) THEN 
RMAX=881 .412109 
GOTO 2 
ENDIF 

3 WRITE (4 , * ) 'The parameters for ' , SFILE, ' were : ' 

WRITE (4, *) ' ’ 

WRITE (4, * ) ' ENERGY= ', ENERGY,’ MeV ’ 

WRITE ( 4 , * ) ' LORENTZ FACTOR= ',G 

WRITE(4,*)' RMS beam angle: ',s,' radians/sqrt (2) . ' 

WR I TE ( 4 , * ) ' IFOIL: ',ifoil 

WRITE ( 4 , *) ' ' 

WRITE ( 4 , * ) ' A value of one indicates that a scattering foil' 
WRITE (4 , * ) ' was present. A value of zero indicates that there' 
WRITE (4 , * ) ' was a target foil only.' 

WRITE (4, *) ’ ' 

WRITE (4 , * ) ' The dielectric constant of the medium was' 

WRITE (4, *) ' ' , XN , ' ^ , XK, ' i . 1 

WRITE ( 4 , *) ' ' 

WRITE ( 4 , * ) ' Angular measurement was over ', ANG, ' radians.' 
WRITE(6, *) 'CHOOSE DESIRED POLARIZATION.' 

WRITE (6, *)’ ENTER APPROPRIATE INTEGER FOR DESIRED POLARIZATION' 
WRITE (6, *) ’ ' 

WRITE(6, *)’l NO POLARIZATION — TOTAL INTENSITY' 

WRITE (6, *)'2 PARALLEL.’ 

WRITE (6, *)'3 PERPENDICULAR.' 

WRITE (6, * ) ' ' 

READ (6,*) POLAR 

IF (ABS (POLAR-1) . LE . IE-6) THEN 

WRITE (4, *) 'TOTAL POLARIZATION WAS SELECTED.' 

ENDIF 

IF (ABS (POLAR-2) . LE. IE-6) THEN 

WRITE ( 4, *) 'PARALLEL POLARIZATION ONLY.' 

ENDIF 
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IF (ABS (POLAR-3) . LE. IE-6) THEN 

WRITE(4, *) 'PERPENDICULAR POLARIZATION ONLY.’ 

END IF 

IF (ABS (IFOIL-1) .LE. IE-6) THEN 
TWOFOIL=0 . 5 

WRITE (6, *) 'SCATTERING FOIL ASSUMED TO BE PRESENT' 

ENDIF 

IF (S *G .GE. 1.0) THEN 

WRITE(6, *) 'WARNING: S MUST BE LESS THAN 1 /GAMMA FOR' 

WRITE (6, *) 'APPROXIMATIONS USED IN I_PERP . ’ 

ENDIF 

* * * *TWFOIL ADJUSTS ASYMMETRIC TERM IN IPAR; WHEN THERE IS A 

* * * * SCATTERING FOIL, SET IFOIL=l, THUS TWFOIL=0.5 

B = SQRT (1 . - (1 . /G) * *2) 

S2= 2.0 * S* *2 
S3=SQRT (PI *S2) 

CS=COS (PSI) 

SN2=1-CS**2 

EPS=XN* *2 - XK* *2 + 2 . *AI *XN*XK 

RPAR= (EPS*CS-CSQRT (EPS-SN2) ) / (EPS*CS+CSQRT (EPS-SN2 ) ) 

RPERP= (CS-CSQRT (EPS-SN2) ) / (CS+CSQRT (EPS-SN2) ) 

RSQ=CABS (RPAR) * *2 
RRL=REAL (RPAR) 

DO 20 1=1, 64 
write (6, * ) ' record= ',i 
DO 10 J=l , 64 

RAD2=REAL ( ( (I-XO) **2) + ( (J-YO) **2) ) 

RAD=SQRT (RAD2 ) 

THETA ( J) =RAD/ SCALER 
IF ( J . LT . 32 ) THEN 

THETA (J) = (-THETA (J) ) 

ENDIF 

IF ( I . GT . 32 . AND . J . GT . 32 ) THEN 
THETA (J)= (-THETA (J) ) 

ENDIF 

X (J) =THETA ( J) 

IF(ABS(I-1) .LT.1E-6) THEN 
X1=X (1) 

ENDIF 

IF (ABS (1-64) . LT . IE- 6 .AND .ABS ( J-64 ) .LT.1E-6) THEN 
X2=X (64) 

ENDIF 
T=THETA (J) 

F = T**2/ (1 . 0/G**2 + T * * 2 ) * *2 
IF (S.EQ. 0.) THEN 

D= (B* *2/PI * *2) * (RSQ*F+ (T/ABS (T) ) ) *RRL*SQRT (F) 

ENDIF 

IF (S . GT . 0 . ) CALL DIST (T, G, B, S, S2, S3, TWFOIL, RPAR, RPERP, DD, D) 
IF (IPRINT .EQ. 1) PRINT *, T, ’ ' , D, ' ',F 

Y ( J) =D 

IF ( S .GT. 0.0) THEN 
Z (J) =DD 

TAU (J) = (DD-D) / (DD t D) 

TAU ( J) =ABS (TAU ( J) ) 
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C* * * “W (I) IS ITOT = 

W (J) =D + DD 
END IF 

IF ( S .EQ. 0.0) THEN 
TAU ( J) =1.0 
Z(J)=0.0 
W (J) =D 
END IF 

IF (Y (J) .GE. YMAX) THEN 
YMAX=Y (J) 

END IF 

IF (Z (J) .GE. PMAX) THEN 
PMAX=Z (J) 

ENDIF 

IF (W (J) .GE.WMAX) THEN 
WMAX=W (J) 

AMAX=ABS (THETA (J) ) 

MAXY=Y (J) 

MAXP=Z (J) 

ENDIF 

IF (ABS (POLAR- 1 ) . LE . IE- 6) THEN 
INT (J) =W ( J) 

ENDIF 

IF (ABS (POLAR-2) . LE. IE-6) THEN 
INT (J) =Y (J) 

ENDIF 

IF (ABS (POLAR-3) . LE . 1 E- 6) THEN 
INT (J) =Z (J) 

ENDIF 

10 CONTINUE 

WRITE (8, REC=I) (INT (COUNT) ,COUNT=l, 64) 

20 CONTINUE 

WRITE (3, 300) ENERGY, S,XDUM, XDUM 
WRITE (3, 300) XI, X2 
WRITE (3, 300) YMAX, WMAX, XDUM 
WRITE<3,*) ROW, COL 

WRITE (4, *) 'MAXIMUM TOTAL INTENSITY WAS ' , WMAX 
WRITE ( 4 , * ) ’ MAX TOTAL INTENSITY OCCURED AT ' , AMAX, ' RADIANS’ 
WRITE(4, *) 'MAX INTENSITY PREDICTED FOR ' , TMAX , * RADIANS' 
WRITE (4, * ) 'MAXIMUM PARALLEL INTENSITY WAS: ' , YMAX 

WRITE(4, *) 'MAXIMUM PERPENDICULAR INTENSITY WAS: ' , PMAX 

WRITE (4, *) 'PARALLEL INTENSITY AT MAXIMUM INTENSITY: ' , MAXY 

WRITE(4, *) 'PERP. INTENSITY AT MAXIMUM INTENSITY: ' , MAXP 

WRITE (6, *) 'NORMALIZING THE INTENSITY VALUES TO 1.0' 

DO 250 NI=1 , 64 

READ (8, REC=NI) ( INT (NW) , NW=1 , 64) 

DO 200 NNW=1 , 64 
IF ( I NORM . EQ . 1 ) THEN 

IF (ABS (POLAR-2) .LE. IE-6) THEN 
INT (NNW) =1 NT (NNW) /YMAX 
GOTO 100 
ENDIF 

IF (ABS (RMAX-0 . ) . LE. 1 .E-6) THEN 
INT (NNW) =INT (NNW) /WMAX 
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100 

c 

c 

c 



300 

200 

222 

250 



ELSE 

INT (NNW) = INT (NNW) /RMAX 
ENDIF 
END1F 

* * * * NORMALIZES IPAR to unity WITH MAX VALUE OF IPAR, YMAX 

* * * *normalizes itot=w(i) to unity with wmax 
'“'normalizes iperp=dd with wmax also 

WRITE (3, 300) X (NNW) , Z (NNW) , Y (NNW) , W (NNW) , TAU (NNW) 

FORMAT (5 (E14.7,1X) ) 

CONTINUE 

WRITE (5, 222) (INT (NWW) , NWW=1, 64) 

FORMAT (64F8 . 4) 

CONTINUE 
CLOSE (2) 

CLOSE (3) 

CLOSE (4) 

CLOSE (5) 

CLOSE (8) 

END 

SUBROUTINE DIST (T, G, B, S, S2, S3, TWFOIL, RPAR, RPERP , IPERP, IPAR) 
T = ANGLE OF OBSERVATION IN RADIANS 
G = LORENTZ FACTOR, GAMMA 
S = RMS SCATTERING ANGLE OR BEAM ANGLE 
D = RESULTING DISTRIBUTION OF INTENSITY 
COMPLEX BETA, Z, W, AI , D1Z , El Z , RPERP , RPAR 
REAL IPERP, IPAR 
DATA PI /3. 141592654/, MO/1/ 

AI = (0 . 0, 1 . 0) 

BETA = 1 . 0/G - AI * T 
Z = BETA/ <1 . 414213562*S) 

IF ( CABS (Z) .LT. 6.1644) THEN 
IF MO = 1, W = ERFC(Z), IF MO = 0, W = ERF (Z) 

CALL CERF (MO, Z , W) 

W=W*CEXP (Z* *2 . ) 

ELSE 



CALL LCERF (Z, W, IER) 

* * * * LCERF RETURNS W=CEXP(Z**2) * ERFC(Z) FOR ABS(Z) . GT . SQRT(38) 

ENDIF 

D1Z = (G - BETA/ S * * 2 . ) * W 
E1Z = (G + BETA/S**2.) * W 
Dl= REAL (D1Z) 

El = REAL (ElZ) 

D2 = (B*S*G) * *2 . * ( (l./S3)**3. * S2*Dl/2. + 1 . / (P I * S3 * * 2 . ) ) 

E2 = B**2. * ( ( 1 . / S 3 ) * * 3 . * S2*El/2. - 1 . / (P I * S3 * * 2 . ) ) 

IPERP = (CABS (RPERP) * *2.) * D2 

IPAR = (CABS (RPAR) * *2 . ) * E2 + TWFOI L* B* * 2 . * REAL (RPAR) * AIMAG (W) 

+ / ( PI* SQRT (PI * S2 ) ) 

RETURN 

END 

SUBROUTINE LCERF (U, W, IER) 

****U IS COMPLEX INPUT VALUE 

****W IS THE VALUE OF CEXP (-U* * 2 ) ERFC (U) OUTPUT 

* * * * REQUIRES ABS ( ARG (U) ) LESS THAN ARG=3*PI/4 

* * * * IF CABS (U) LESS THAN SQRT (38) THIS ROUTINE IS NOT VALID 
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c 

c 

c 



,... 1ER RETURNS AS 999 IF CONDITION ON ^10)^ IS VIOLATED 
1 1 * * ir^THE^SUBROUTINE^ CONDITIONS ARE VIOLATED RETURNS W-«IER..ER. 

DAT^pI/^^^ 

IER=0 

X=REAL(U) 

Y=AIMAG (U) 

TARG=ATAN2 (Y, X) 

IF (ABS (TARG) • GT . ARo) TH 
W= (999 , 999) 

IER=999 

RETURN 

iff CABS <UI -LT. 6.1644, THEN 
W= (-999 . » -999 . ) 
lER=-999 
return 

END IF 

US=2 . 0 *U*U 
USUM= (1 . 0, 0 .0) 

UPROD= (1 . 0 , 0 . 0) 

D ° 10 UPROD=-UPROD‘(2* 1-1) /US 

USUM=USUPU UPROD 
1f) CONTINUE 

W=USUM/ (U*RTPI) 

RETURN 

SUBROUTINE CERF (MO ; _A, ^ . “’**"*" 

TORT^* SOB^TINB^rOR .COMPLEX ^ ERROR ^ FONCT IO« 

« nunber, 

5 : FUNCTION (COMPLEX NUMBER, 

COMPLEX A, W C7<?} EF (2) 

DIMENSION AZ(2), Q2(2), ' QF( 2), CD(18), CE ( 1 ® \ n / 

DMA S(9)a‘J 

iillliilllliilii! ' 

nRTA CE(1) /8.15723083324096E 2/,C ( , x _ 3321 967 083624 5E-1 / 

DMA CE (5) / I 1569039287 8 957E- 1/' CE ( 6) /9 ^ ^04 538609872E-2/ 

srMU;:S;:«/L3»«6 SM 7,2» B -2/ 
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DATA CE (11) /2. 0554809947019 3E -2/, CE (12) /I . 07847 403887506E-2/ 
DATA CE (13) /4 .55 6348922 1421 9E-3/,CE( 14) / 1.4 3984 4 58138 925E-3/ 
DATA CE (15) /3. 07056139834171 E- 4/, CE (16) /3 . 78 15654 1 1 6854 IE-5/ 
DATA CE (17) /2.05173509616121E-6/, CE ( 18 ) /2 . 63564823682747E-8/ 



AZ ( 1 ) =REAL (A) 

AZ (2) =AIMAG (A) 

001 ZS=AZ (1) *AZ (1) +AZ (2) *AZ (2) 

SZ (1) =AZ (1) *AZ (1) -AZ (2) *AZ (2) 

SZ (2) =2 . 0*AZ (1) *AZ (2) 

QZ ( 1 ) =+AZ ( 1 ) 

QZ (2) =+AZ (2) 

SN=+ 1 . 0 

1 F ( AZ ( 1 ) ) 002, 003, 003 

002 QZ ( 1 ) =-AZ ( 1 ) 

QZ (2) =-AZ (2) 

SN=-1.0 

003 IF (ZS-1 . 0) 014, 014, 004 

004 IF (ZS-38 . 0) 005, 006, 006 

005 IF (SZ (1) -tO. 064 *SZ (2) *5Z(2))014,014,012 

ASYMPTOTIC EXPANSION 

006 SM ( 1 ) =0 . 0 
SM (2 ) =0 . 0 

QF (1) =+SZ (1) / (ZS*ZS) 

QF (2) =-SZ (2) / <ZS*ZS) 

QM=0 .5 64 18 958 354 7 75 6* EXP (-SZ(l) ) 

TS (1) =QM*COS (-SZ (2) ) 

TS (2) =QM*SIN (-SZ (2) ) 

TM(1) =+ (TS (1) *QZ(1) ITS (2) *QZ (2) ) /ZS 
TM(2) =- (TS (1) *QZ (2) -TS (2) *QZ (1 ) ) / ZS 
PM=-0 . 5 
GO TO 008 

007 PM=PM+ 1 . 0 

TS ( 1 ) =TM ( 1 ) *QF ( 1 ) -TM (2 ) *QF(2) 

TS (2) =TM ( 1 ) *QF (2) i-TM (2) *QF (1) 

TM ( 1 ) =-PM*TS (1) 

TM (2) =-PM*TS (2) 

IF (ABS (SM ( 1 ) ) +ABS (TM ( 1 ) ) . NE . ABS (SM (1 ) ) ) GO TO 008 
IF (ABS (SM (2) ) +ABS (TM (2) ) .EQ.ABS (SM (2) ) ) GO TO 009 

008 SM ( 1 ) =SM ( 1 ) +TM ( 1 ) 

SM (2) =SM (2)4 TM (2) 

1F(PM.LT.36.5) GO TO 007 

009 IF (QZ(1) .LT.0.01) GO TO 019 
EF ( 1 ) =SM ( 1 ) 

EF (2) =SM (2) 

010 IF (MO . EQ . 0 ) GO TO Oil 
W=CMPLX (EF (1) , EF (2) ) 

IF (SN . EQ .1.0) RETURN 
W=CMPLX (2 . 0-EF (1) , -EF (2) ) 

RETURN 

011 EF (1) =SN*(1 • 0-EF (1) ) 

EF (2) =-SN*EF (2) 
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W=CMP LX (EF ( 1 ) , EF (2 ) ) 
RETURN 



RATIONAL FUNCTION APPROXIMATION 

012 SM ( 1 ) =0 . 0 
SM (2 ) =0 . 0 

QM=0 .564189583547756’*EXP (-SZ(I) ) 

TS (1) =QM*COS (-SZ (2) ) 

TS (2) =QM*SIN (-SZ (2) ) 

QF (1) =TS (1) *QZ (1) -TS (2) *QZ (2) 

QF (2) =TS (1 ) *QZ (2) +TS (2) *QZ(1) 

DO 013 1=1, 18 
TS (1) =SZ (1) +CD (I) 

TS (2) =SZ (2) 

SS=TS (1) *TS (1) +TS (2) *TS (2) 

TM ( 1 ) =+CE (I) *TS (1) /SS 
TM ( 2 ) =-CE (I) *TS (2) /SS 
SM ( 1 ) =SM ( 1 ) +TM ( 1 ) 

SM (2) =SM (2) +TM (2 ) 

013 CONTINUE 

EF ( 1 ) =QF ( 1 ) * SM ( 1 ) -QF (2 ) *SM(2) 

EF ( 2 ) =QF ( 1 ) * SM (2) +QF (2) *SM ( 1 ) 

GO TO 010 

TAYLOR SERIES EXPANSION 

01-1 TM(1) =1 . 1 2837 91 6709551 3*QZ (1) 

TM (2) =1 .1283791 67095513 *QZ (2) 

SM ( 1 ) =TM ( 1 ) 

SM (2) =TM (2) 

PM=0 . 0 

015 PM=PM+ 1 . 0 
DM=2 . 0*PM+1 . 0 

TS (1) =TM (1) *SZ (1) -TM (2) *SZ (2) 

TS (2 ) =TM ( 1 ) * SZ (2) 4 TM ( 2 ) *SZ(1) 

TM (1) = — TS (1) /PM 
TM ( 2 ) =-TS (2) /PM 
TS (1) =TM (1) /DM 
TS (2) =TM (2) /DM 

IF (ABS (SM ( 1 ) ) -t ABS (TS ( 1 ) ) . NE . ABS (SM (1 ) ) ) GO TO 016 
IF (ABS (SM ( 2) ) +ABS (TS (2 ) ) . EQ . ABS (SM (2 ) ) ) GO TO 017 

016 SM ( 1 ) =SM ( 1 ) +TS (1) 

SM (2 ) =SM (2 ) 4TS(2) 

GO TO 015 

017 IF (MO.NE.O)GO TO 018 
W=CMPLX (SN* SM ( 1 ) , SN*SM(2) ) 

RETURN 

018 EF (1 ) =1 . 0-SN* SM ( 1 ) 

EF (2) =-SN* SM ( 2 ) 

W=CMPLX (EF (1) , EF (2) ) 

RETURN 

MODIFIED ASYMPTOTIC EXPANSION 
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019 SN=-SN 

GO TO 017 

r>in 



B. TWO-FOIL INTERFEROMETER PROGRAM 



PROGRAM COHER3 

C REVISION OF TEM2 , COHER1 BY D.W. RULE 
C REVISED BY WIL LONGSTAFF, NPGS AUG, 1988 
C PLOTS NORMALIZED TR DATA IN THREE DIMENSIONS 

C INCLUDES CHANGES MADE BY D.W. RULE FOR COHERENT CONTRIBUTION 

C OF A CLEAR FRONT FOIL 

C* ' * * * D. W. RULE, NSWC MAR. 1987 REVISION OF TINTF 

C* ' *CALC . INTENSITY DIST. AS FUNC. OF ANGLE OF TRANSITION RAD. 

C CALCS. UNPOLARIZED DIST. AND THE PARALLEL & PERP COMPS. 

C ASSUMING GAUSSIAN BEAM SCATTERING ANGLE FOR 
C WARTSKI 1 S INTERFEROMETER 
C 

C G = LORENTZ FACTOR, S=SX=PROJECTED RMS BEAM ANGLE IN 
C RADIANS IN X-Z PLANE 

C SY=PROJECTED RMS BEAM ANGLE IN Y-Z PLANE IN RADIANS 
C T = ANGLE OF OBSERVATION IN RADIANS MEAS . WRT 45 DEGS . 

C INCLUDES BANDWIDTH OF WAVELENGTHS WLMIN TO WLMAX 
CHARACTER DATFI L*50, CH* 1 , SFILE* 15 , DFILE* 1 5 
COMPLEX EPS, AI, 0, EPSIN, EPCOS , ZF, RTQ, Z2, Z4 , RPARF , RPERF 
COMPLEX RPARB, RPERB, EPSF, TCOEF, SDELT, RPARFSQ 
INTEGER AC, BC, NNW, ROW, COL, XO, YO, COUNT, WC, LW, NC, NW, NI, NWW 
INTEGER POLAR, CHOICE, ENTRY, N, OPTION 

REAL THETA ( 64 ) , RAD2, RAD, SCALER, ZMIN, ZMAX, ZCEN, INT(64) 

REAL X (64) , TAU (64) , I PERP (64) , IPAR<64) , W ( 64 ) , IPP , I PR 
REAL IPPF, IPPB, IPRF, IPRB, KOEFF, C ( 4 , 64 ) , DX ( 64 ) , DY ( 64 ) 

REAL TMAX, AMAX, ANG, MAXY, MAXP , NMBR, PMAX, THMX, RMAX, QMAX 
C * * * ENERGY = BEAM ENERGY IN MEV 

C * * *XNB, XKB ARE REAL AND IMAG OPTICAL CONSTANTS OF METALIC (B)ACK FO 
C * * * XNF , XKF ARE OPTICAL CONSTANTS OF TRANSPARENT (F) RONT FOIL 
C * * *XL IS SEPARATION OF INTERFEROMETER IN CM 
C ***PSI IS ANGLE OF BEAM WRT NORMAL TO THE FOIL, PSI=PI/4 

DATA PI/3. 141592654/, PSI/0 . 7853982/, THMX /0 . /, MAXY/0 . /MAXP/0 . / 

DATA IFILE/3/ , I PRINT/0/, INORM/1/, IFOIL/O/, SX/ .000210/ 

DATA SY/ . 000120/, XL/1 . 2/, DFOIL/ . 005/WLMIN/4500 . / , DELT/7 . 79E-4/ 

DATA WLMAX/ 4 5 00 . 0/, XNF/1 . 4 8/, XKF/0 . / , XNB/ . 618/, XKB/5 .4 7/ 

DATA KK/3/,FFl/l .0/, TWFOIL/1 .0/, XO/32/, YO/32/, ROW/64/ 

DATA COL/64/, SCALER/620. 25841/, ZMIN/0.0/, ZMAX/1 .0/, ZCEN/0./ 

DATA RMAX / 0 . 0 / , QMAX / 0 . 0 / , PMAX / 0 . 0 / , INT/64*0.0/ 

WRITE(6,*)’ ENTER NAME OF GRAPHICS DATA FILE.’ 

WRITE (6,*)' EXAMPLE: GRAPH . DAT ' 

WRITE (6, *) ’ ' 

READ (6, 7) SFILE 

WRITE(6,*)’ ENTER NAME OF CORRESPONDING PARAMETER FILE.' 

WRITE (6,*)’ EXAMPLE: PRMTRS.DAT' 

WRITE (6, *) ’ ’ 

READ (6, 7) DFILE 
7 FORMAT (A15) 

OPEN ( 2, F I LE=' SCRATCH' , STATUS= ' UNKNOWN ’ , ACCESS= ' DIRECT ' , 

+ RECL=256) 

OPEN (3, FILE= ' TOUT . DAT ' , STATUS= ' UNKNOWN ’ ) 

OPEN (4, FILE=SFILE, FORM= ’ FORMATTED ' ) 

OPEN (5, FILE=DFILE, STATUS= ' UNKNOWN ' ) 

WRITE (4 , 123) ROW, COL, ZMIN, ZMAX, ZCEN 
122 FORMAT (2 1 3 , 3F8 . 4 ) 
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WRITE (6, *) ’WHAT IS THE DESIRED BEAM ENERGY IN MEV? 
WRITE (6, *) 'PLEASE INCLUDE A DECIMAL PLACE.' 

WRITE (6, *) 1 ' 

READ ( 6 , * ) ENERGY 
G=1 . 0 + ENERGY/0.511 
THMX= 1 . /G 



ANG =31 
WRITE (6, *) ’The 
WRITE ( 6, *) 1 ' 

WRITE (6, 

WRITE (6, 

WRITE (6, 

WRITE ( 6, 

WRITE (6, 

WRITE (6, 

WRITE (6, 

WRITE (6, 

WRITE (6, 

WRITE ( 6, 



/SCALER 

current 



parameters are: 



*) 

) 

) 

') 
* ) 
*) 



angle in 
angle in 
' , IFOI L, ’ 



the 
the 
. A 



plane 
plane 
value of 1 



X-Z 

Y-Z 



is 

is 



1 1 RMS beam 
' RMS beam 
'2 IFOI L IS 
' that 
' zero 

'3 The optical constants 
’ are 1 , XNB, '+ ' , XKB, ' i . ' 

'4 The optical constants 
’ are ' , XNF, ’ + ’ , XKF, ’ i . ' 

'5 Angular measurement is over ' , ANG, ' radians 
IF (ABS (RMAX-0. ) . LE.1E-6) THEN 

WRITE (6,*)' 6 Normalization is to the maximum total' 



radians 

radians 



\SX, 

*,SY, 

indicates ' 

a scattering foil is present. A value of' 
indicates absence of scattering foil.' 

for the back metalic foil 

for the front clear foil* 



WRITE ( 6, A )' intensity for 
WRITE (6, *) ’The option to 
'intensity for 
'In all casis, 
'normalized to 
'energy entered 



)' Normalization will be to 100 MeV.' 

Minimum wavelength is *,WLMIN,' angstroms. 
Maximum wavelength is ',WLMAX,' angstroms. 



Foil thickness is ',DFOIL, ' 
Distance between foils is 1 
ACCEPT CURRENT PARAMETERS' 



WRITE ( 6 , 

WRITE (6, 

WRITE (6, 

WRITE ( 6, 

ELSE 

WRITE (6, 

END IF 

WRITE (6, *) '7 
WRITE (6, *) ' 

WRITE (6, *) '8 
WRITE(6, *) ’9 
WRITE (6, * ) ’ 10 
WRITE ( 6, *) ' ' 

READ ( 6 , * ) CHOICE 
IF (ABS (CHOICE-10) .LE.1E-6) THEN 
GOTO 3 
END IF 

IF (ABS (CHOICE-1) .LE.1E-6) THEN 
WRITE (6, *) 'ENTER NEW VALUE FOR 
WRITE ( 6, *) ' ’ 

READ ( 6, * ) SX 
WRITE (6, * ) 'ENTER 
WRITE ( 6 , *) ’ ' 

READ ( 6 , * ) SY 
GOTO 2 
END IF 

IF (ABS (CHOICE-2) .LE.1E 
WRITE (6, * ) 'ENTER NEW 
WRITE ( 6 , *) ' ' 

READ ( 6 , * ) IFOI L 
GOTO 2 



the energy entered for this run. 
normalize to the maximum total' 
an energy of 100 MeV. ' 
parallel intensity will be' 
its own maximum value for the' 
for this run. ' 



centimeters . 

XL,' centemeters. 



X-Z PLANE RMS BEAM ANGLE’ 



NEW VALUE FOR Y-Z PLANE RMS BEAM ANGLE’ 



■6) THEN 
VALUE FOR 



I FOIL . 
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ENDIF 

IF (ABS (CHOICE-3) . LE . 1 E- 6 ) THEN 

WRITE ( 6, * ) ' ENTER THE REAL VALUE FOR THE BACK FOIL CONSTANT 
WRITE (6, *) ’ ' 

READ (6,*) XNB 

WRITE(6, *) 'ENTER THE IMAGINARY VALUE FOR THE BACK FOIL’ 
WRITE (6, *) ’CONSTANT. DO NOT INCLUDE THE i.' 

WRITE (6, *) ’ ’ 

READ (6,*) XKB 
GOTO 2 
ENDIF 

IF (ABS (CHOICE- 4 ) . LE. IE- 6) THEN 

WRITE (6, *) ’ENTER THE REAL VALUE FOR THE FRONT FOIL CONSTANT 
WRITE (6, *) ’ ' 

READ (6,*) XNF 

WRITE (6, *) 'ENTER THE IMAGINARY VALUE FOR THE FRONT FOIL’ 
WRITE (6, *) 'CONSTANT. DO NOT INCLUDE THE i.' 

WRITE (6, *) ' ' 

READ (6, * ) XKF 
GOTO 2 
ENDIF 

IF (ABS (CHOICE-5) .LE.1E-6) THEN 
WRITE(6,*) 'Choose option.’ 

WRITE(6,*)’l Angle that is N^ ( 1/Lorentz factor)’ 
WRITE(6,*)'2 Enter angle in radians manually.' 

WRITE (6, * ) ’ ' 

READ (6,*) OPTION 

IF (ABS (OPTION-1) .LE. IE-6) THEN 

WRITE ( 6, *)' Enter N, for N* (1/Lorentz factor)’ 

WRITE (6, *) ' ’ 

READ (6,*) NMBR 
ANG=NMBR/G 

SCALER=3 1 . / (SIN (ANG) ) 

ENDIF 

IF (ABS (OPTION-2) .LE.1E-6) THEN 

WRI TE ( 6, * ) ' Enter desired angle in radians.' 

WRITE (6, *) ' ' 

READ (6,*) ANG 
SCALER=3 1 . / (SIN (ANG) ) 

ENDIF 
GOTO 2 
ENDIF 

IF (ABS (CHOICE-6) .LE.1E-6) THEN 
RMAX=881 .412109 
GOTO 2 
ENDIF 

IF (ABS (CHOICE-7) .LE. IE-6) THEN 

WRITE(6,*) 'ENTER NEW MINIMUM WAVELENGTH. ' 

WRITE (6, *) ’ ’ 

READ (6,*) WLMIN 

WRITE (6, *)’ ENTER NEW MAXIMUM WAVELENGTH. ' 

WRITE (6, *) ’ ’ 

READ ( 6 , * ) WLMAX 
GOTO 2 
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ENDIF 

IF (ABS (CHOICE-8) .LE. IE-6) THEN 

WRITE (6, * ) ' ENTER VALUE FOR FOIL THICKNESS IN CENTIMETERS.' 
WRITE (6, *) ' ’ 

READ (6,*) DFOIL 
GOTO 2 
ENDIF 

IF (ABS (CHOICE-9) .LE. IE-6) THEN 

WRITE (6, *) 'ENTER DISTANCE BETWEEN FOILS IN CENTIMETERS.’ 
WRITE (6, *) ’ ’ 

READ (6,*) XL 
GOTO 2 
ENDIF 

3 WRITE (5, *)' The parameters for • , SFILE, ’ were:' 

WRITE (5, *) ' ' 

WRITE (5, *) 'ENERGY= ENERGY,’ MeV’ 

WRITE (5, * ) ' LORENTZ FACTOR= \G 

WRITE ( 5 , * ) ’ RMS beam angle in the X-Z plane: ' , SX 
WRITE (5,*)’ RMS beam angle in the Y-Z plane: ' , SY 
WRITE (5, * ) ’ IFOIL: ’.IFOIL,’ A VALUE OF 1 INDICATES’ 

WRITE (5 , * ) ' THAT A SCATTERING FOIL WAS PRESENT. A VALUE OF' 
WRITE (5, *)' ZERO INDICATES ABSENCE OF SCATTERING FOIL.’ 

WRITE (5,*) 'THE DIELECTRIC CONSTANTS FOR THE BACK METALIC FOIL 
WRITE (5, * ) 1 ' , XNB, ' , XKB, ' i . ’ 

WRITE(5,*) 'THE DIELECTRIC CONSTANTS FOR THE FRONT CLEAR FOIL 
WRITE<5,*)' ' , XNF, ’ + ' , XKF, ’ i . ' 

WRITE (5, * ) ' Angular measurement was over ',ANG, ' radians.' 

IF (ABS (RMAX-0 . ) .LE.1E-6) THEN 

WRITE(5,*)' Normalization to energy selected.’ 

ELSE 

WRITE (5, *)' Intensity was normalized to 100 MeV.’ 

ENDIF 

WRITE (5, *) ’ MINIMUM WAVELENGTH: \WLMIN 

WRITE(5,*)' MAXIMUM WAVELENGTH: ' , WLMAX 

WFlTE(5,*)'Foil thickness: ', DFOI L, ' centimeters . ' 

WRITE (5,*)' Distance between foils: ' , XL, ’ centimeters.' 

WRITE (6, *) 'CHOOSE THE DESIRED POLARIZATION.’ 

WRITE(6, *) 'ENTER APPROPRIATE INTEGER FOR DESIRED POLARIZATION 
WRITE (6, * ) ' ’ 

WRITE(6,*)'l NO POLORI ZATION--TOTAL INTENSITY.' 

WRITE (6,*)' 2 PARALLEL.’ 

WRITE (6,*) '3 PERPENDICULAR.’ 

WRITE (6, *) ’ ’ 

READ (6,*) POLAR 

IF (ABS (POLAR-1) . LE . ] E-6 ) THEN 
WRITE(5, *) 'TOTAL POLARIZATION WAS SELECTED' 

ENDIF 

IF (ABS (POLAR-2) . LE. IE-6) THEN 
WRITE (5, *)' PARALLEL POLARIZATION ONLY’ 

ENDIF 

IF (ABS (POLAR-3) .LE. IE-6) THEN 
WRITE (5, *)' PERPENDICULAR POLARIZATION ONLY' 

ENDIF 

c****ffl is factor for incoh. addn . of single foil intensities 
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****from clear foil ffl=1.5 for looking thru foil 
CALCULATING THE PROGRAM CONSTANTS 
S=SX 

VDC = SORT (1.0 - ( 1 . 0/G) **2.0 ) 

AI= (0 . 0, 1.0) 

IF ( IFOIL .EQ. 1) THEN 
TWFOIL=0 . 5 
ENDIF 

BMIN=P I * XL/ (2 . 0*VDC*WLMAX* 1 .0E-08) 

BMAX=BMIN * WLMAX /WLMIN 
IF ( WLMAX .EQ. WLMIN ) THEN 
BMAX=BMIN 
ENDIF 

S2=2 . 0*S* *2 
S3=SQRT (PI *S2) 

DELB = BMAX- BMIN 

IF ( S .GT. 0.0 .AND. BMAX . EQ . BMIN) THEN 
DELB =1 .0 
B=BMIN 

Q= 1 . 0/S2-AI *2 . 0*B 
RTQ=CSQRT (Q) 

Z2=AI *RTQ/G 
Z4=Q/G 

END IF 

CALL REFLEC (PSI, XNF, XKF, RPARFS, RPERFS, RPARF, RPERF) 

CALL REFLEC (PSI, XNB, XKB, RPARBS, RPERBS, RPARB, RPERB) 

***RPARFS, RPERFS ARE PARALLEL & PERP. REFLECTIVITIES OF FRONT FOIL 

* * * WHICH IS TRANSPARENT, RPARBS AND RPERBS ARE FOR BACK METALIC FOIL 

* * * RPARF RPERF AND RPARB RPERB ARE THE CORRESPONDING REFLECTION COEFFS 

KOEFF = VDC* *2.0 *RPERBS*SY **2.0/PI**2.0 
IF (S .EQ. 0.0) THEN 
B=BMI N 
ENDIF 

DO 20 1 = 1, 64 

WRITE (6, *) 'RECORD= ',1 
DO 10 J=l, 64 

RAD2=REAL ( ( (I-XO) **2) + ( (J-YO) **2) ) 

RAD=SQRT (RAD2) 

THETA < J) =RAD/ SCALER 
I F ( J . LT . 32 ) THEN 

THETA ( J) = (-THETA ( J) ) 

ENDIF 

X (J) =THETA (J) 

IF (ABS ( I- 1 ) . LT . IE- 6 . AND . ABS ( J - 1 ) .LT.lE-6) THEN 
TMIN=X (1) 

TMAX= (-TMIN) 

WRITE (IFILE, 15) ENERGY, G, SX, SY, XL 
WRITE (IFILE, 15) TMIN, TMAX, DELT, WLMAX 
WRITE (IFILE, 15) XNF, XKF, XNB, XKB, WLMIN 
WRITE (IFILE, *) ROW, COL 
15 FORMAT (5 (E14 .7, IX) ) 

ENDIF 
T=THETA ( J) 

IF (S .GT. 0.0) THEN 
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C * *CALC . SINGLE FOIL PATTERNS IPPF (PERP) f. 1PRF (PAR) FROM (F) RONT FOIL 

CALL DIST (T, G, VDC, SX, SY, S2, S3, TWFOIL, RPARFS, RPERFS, 

-tRPARF, IPPF, IPRF) 

C * * *CALC . SINGLE FOIL PATTERNS IPPB (PERP) & IPRB (PAR) FROM (B) ACK FOIL 

TFL =1.0 

CALL DIST (T, G, VDC, SX, SY, S2, S3, TFL, RPARBS, RPERBS , RPARB, 

+ IPPB , IPRB) 

IF (BMAX .GT. BM1N) THEN 

CALL BANDIN (KK, T, G, S, S2, BMIN, BMAX, DPAR, DPERF) 

ELSE 

CALL DIST2 (T,G, S, Q, B, S2, RTQ, Z2, Z1 , DPAR, DPERF) 

END IF 

IPAR(J) = 2.0* (ffl*IPRB- RPARBS* DPAR/ (DELB *PI**2.0) )+2.*lPRF 

IPERP(J) = 2.0* ( f f 1 * IPPB-KOEFF*DPERF/DELB ) +2.*IPPF 

* *NOTE FACTOR OF 1.5 IN IPAR & IPERP TO SIMULATE 1ST SURFACE 
“CONTRIBUTION WHICH IS INCOHERENT 

IPAR, IPERP, OR W IS INTENSITY/ (UNIT FREQ . , UNIT SOLID ANGLE) 

* * * IN UNITS OF (CHARGE* *2) / (VEL OF LIGHT) 

ELSE 

FT=1 . 0/G* *2 + T* *2 

D3~ (T/FT) * * 2 

IPAR ( J) = / l . 0 * RPARBS / P 1**2 * D3 *SIN (B*FT) * * 2 

IPERP ( J) -A .0 *KOEFF * ( SIN(B*FT) / FT )**2 

END IF 

W(J) = IPAR (J) * IPERP (J) 

C * * * I F S=0.0, D IS CALCULATED FROM WARTSKI’S EXPRESSION 

IF ( IPRINT .EO. 1 ) PRINT *, T,' ',D 

IF (ABS ( POLAR- 1 ) . LE . 1E~6 ) THEN 
INT (J) =W (J) 

END IF 

IF (ABS (POLAR-2) . LE . IE-6) THEN 
INT ( J) =IPAR (J) 

ENDIF 

IF (ABS (POLAR-3) .LE. IE-6) THEN 
INT ( J) = I PERP (J) 

ENDIF 

10 CONTINUE 

WRITE (2, REC=I ) (INT (COUNT) , COUNT=l , 61) 

CALL NORM (T, IPAR, IPERP, W, TAU, I NORM, RMAX, QMAX , PMAX , AMAX , MAXY , 

' MAXP) 

DO 550 BC=1 , 61 

WRITE (IFILE, 33) X(BC) , IPERP(BC) , IPAR(BC) ,W(BC) ,TAU(BC) 

33 FORMAT (5 (Ell . 7, IX) ) 

550 CONTINUE 

20 CONTINUE 

DO 11 AC=1, 61 

READ (2, REC=AC) ( INT (WC) , WC=1 , 61 ) 

C * * * * NORM REPLACES w WITH IPAR + IPERP, I.E. THE UNPOLARIZED INT 

IF ( I NORM . EQ. 1) THEN 
DO 200 NW=1 , 61 

IF (ABS (POLAR-2) . LE. IE-6) THEN 
INT (NW) =INT (NW) /OMAX 
GOTO 200 
ENDIF 
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200 



INT (NW) = I NT (NW) /RMAX 
CONTINUE 

ENDIF 



C 



29 



4 1 



WRITE (A , 29) (INT (NW) , NW=1, 64) 
FORMAT ( 64 F 8. 4) 

WRITE (6, * ) (INT (NW) , NW=I, 64) 
CONTINUE 



WRITE (5, *) 'MAXIMUM TOTAL INTENSITY WAS ' , RMAX 

WRITE (5, *) 'MAX TOTAL INTENSITY OCCURED AT AMAX, ' RADIANS.' 
WRITE (5, *) 'MAX INTENSITY PREDICTED FOR ' , THMX, ' RADIANS.’ 
WRITE (5, *) 'MAXIMUM PARALLEL INTENSITY WAS: ’ , QMAX 

WRITE (5,*) 'MAXIMUM PERPENDICULAR INTENSITY WAS: ’ , PMAX 

WRITE (5 , * ) ' PARALLEL INTENSITY AT MAXIMUM INTENSITY: ’ , MAXY 

WRITE (5, * ) ’ PERP INTENSITY AT MAXIMUM INTENSITY: ’ , MAXP 



END 

SUBROUTINE DIST2 (T, G, S, Q, B, S2 , RTQ, 7,2, ZA, DPAR, DPERP) 

T=ANGLE OF OBSERVATION, THETA RADIANS 
G= LORENTZ FACTOR 

S=SX= RMS BEAM DIVERGENCE ANGLE PROJECTED IN THE X-Z PLANE 
Q=COMPLEX FUNCTION OF S, INTERFEROMETER SPACING L, AND WAVELENGTH 
DPAR AND DPERP ARE THE INTERFERENCE TERMS IN THE DISTRIBUTION 
COMPLEX Q, AI, WC, WD, ZC, ZD, Zl, Z2, Z3, Z4, ZIQC, ZIQD, ZEX, 

4 ZD2, RTQ, EZEX, ZPERP 

DATA PI/3. 141592654/, MO/1/ 

AI= (0 . 0, 1.0) 

Zl=-T/ (S2*RTQ) 

ZC=-AI * (Z1+Z2) 

ZD=AI * (Z1-Z2) 

Z3=-AI*T/S2 
Z5=G-2 . 0*Z4 
ZIQC=Z3-Z4 
ZIQD=Z34 Z4 

ZEX=AI * 2 . 0 *B* (1 . 0/G**2 + T**2) -4 .0* (B*T) * *2/Q 
EZEX=CEXP (ZEX) 

CALL WERFC (ZC, ZD, WC, WD) 

ZD2=0.5*SQRT (PI/S2) *EZEX* ( (G/2. + ZIQD) * WD 

+ +(G/2. - ZIQC) * WC - 2. *RTQ/SQRT (PI) ) 

DPAR = REAL ( ZD2 ) 

ZPERP=0 .25* SQRT (PI/S2) *G**2*EZEX* 

4 ( WD* (G-2 . 0 *RTQ*ZD) 4 WC* (G-2 . 0*RTQ*ZC) 4 4 . 0*RTQ/SQRT (PI) ) 

DPERP = REAL (ZPERP) 

RETURN 

END 

SUBROUTINE WERFC ( ZC , ZD, WC , WD ) 

COMPLEX ZC , ZD, WC, WD, U, W 

MO=l 

IER=0 

DO 10 1=1,2 



IF (I .EQ. 1) THEN 
U=ZC 

ELSE 

U=ZD 

ENDIF 



STOP 
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IF ( CABS (U ) .LT. 6.1644) THEN 
CALL CERF (MO, U, W) 

C * * * * CERF RETURNS W=ERFC(U) IF MO=l, SEE DAHLGREN MATHLIB TR84-14 3 

W=W* CEXP (U * * 2 . ) 

ELSE 

C * * * * LCERF RETURNS W=CEXP (U* * 2 . ) ERFC (U) FOR ABS(U) .GT. SQRT(38) 

CALL LCERF (U,W, IER) 

IF ( I ER .NE. 0) PRINT *,'ARG(U) OUT OF BOUNDS *, 

+ ’IN CALL TO LCERF’ 

ENDIF 

IF (I .EQ. 1) THEN 
WC=W 

ELSE 

WD=W 

ENDIF 

10 CONTINUE 

RETURN 
END 

SUBROUTINE LCERF (U, W, IER) 

**“U IS COMPLEX INPUT VALUE 

* * * * W IS THE VALUE OF CEXP (U* * 2) ERFC (U) OUTPUT 

* * * * REQUIRES ABS ( ARG(U) ) LESS THAN ARG=3*PI/4 

* * * * IF CABS (U) LESS THAN SORT (38) THIS ROUTINE IS NOT VALID 

* * * * IER RETURNS AS 999 IF CONDITION ON ARG(U) IS VIOLATED 

* * * * IER RETURNED AS -999 IF CABS (U) IS TOO SMALL 

* * * * I F THE SUBROUTINE CONDITIONS ARE VIOLATED RETURNS W=(IER, IER) 
COMPLEX U, W, US, USUM, UPROD 

DATA PI/3. 141592654/, ARG/2. 3562/, RTPI/1 .7724 538509055/ 

IER=0 
X=REAL (U) 

Y=AIMAG (U) 

TARG=ATAN2 (Y, X) 

IF (ABS(TARG) . GT . ARG) THEN 
W= (999, 999) 

IER=999 
RETURN 

ENDIF 

IF ( CABS (U) .LT. 6.1644) THEN 
W= (-999 . , -999 . ) 

IER=- 999 
RETURN 

ENDIF 

US=2 . 0*U*U 
USUM=(1. 0,0.0) 

UPROD= (1 .0,0.0) 

DO 10 1=1, 37 

UPROD=-UPROD* (2*1-1) /US 
USUM=USUM+ UPROD 
10 CONTINUE 

W=USUM/ (U* RTP I ) 

RETURN 
END 

SUBROUTINE REFLEC (PSI, XN, XK, RPARS, RPERS, RPAR, RPER) 

COMPLEX AI, EPS, RPAR, RPER 
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AI= (0. 0, 1 . 0) 

CS=COS (PSI) 

SN2=1-CS * * 2 . 

EPS=XN* *2 . - XK* * 2 . + 2 . * AI * XN*XK 

RPAR= (EPS *CS~CSQRT (EPS-SN2) ) / (EPS*CS+CSQRT (EPS-SN2 ) ) 

RPER= (CS-CSQRT (EPS-SN2) ) / (CS+CSQRT (EPS-SN2) ) 

RPARS=CABS (RPAR) * *2 . 

RPERS=CABS (RPER) * *2 . 

RETURN 

END 

SUBROUTINE DIST (T, G, B, SX, SY, S2, S3, TWFOIL, RPARS , RPERS , 

+ RPAR, IPERP, IPAR) 

C T = ANGLE OF OBSERVATION IN RADIANS 

C G = LORENTZ FACTOR, GAMMA 

C S=SX=RMS SCATTERING ANGLE OR BEAM ANGLE PROJECTED IN X-Z PLANE 

C SY=RMS SCATTERING ANGLE OR BEAM ANGLE PROJECTED IN Y-Z PLANE 

C IPERP AND IPAR ARE RESULTING DISTRIBUTION OF INTENSITIES 

C B = VDC IN MAIN PROGRAM, VDC = V/C, I.E. VEL. 

C I OVER SPEED OF LIGHT 

COMPLEX BETA, Z , W, A I , D 1 Z , El Z , RPERP, RPAR 

REAL IPERP, IPAR 

DATA PI /3. 141592654/, MO/1/ 

AI = (0.0, 1 .0) 

S=SX 

BETA = 1 . 0/G - AI * T 
Z = BETA/ (1 .414213562*S) 

IF ( CABS (Z) .LT. 6.1644) THEN 
C IF MO = 1, W = ERFC(Z), IF MO = 0, W = ERF (Z) 

CALL CERF (MO, Z, W) 

W=W*CEXP (Z* * 2 . ) 

ELSE 

CALL LCERF ( Z, W, IER) 

C * * * * LCERF ETURNS W=CEXP(Z**2) * ERFC(Z) FOR ABS (Z) . GT . SQRT(38) 

END IF 

D1Z = (G - BETA/ S * * 2 . ) * W 
E1Z = (G + BETA/S * * 2 . ) * W 
Dl = REAL (D1 Z) 

El = REAL (E1Z) 

D2 = (B * S Y *G) * * 2 . * ( (l./S3)**3. * S2*Dl/2. + 

+ 1 . / (PI *S3* *2 . ) ) 

E2 = B * * 2 . * ( (l./S3)**3. * S2*El/2 . - 1 . / (PI 'S3* *2 . ) ) 

IPERP = RPERS * D2 /TWFOIL 

IPAR = RPARS * E2 / TWFOIL + B* * 2 . * REAL (RPAR) * A IMAG (W) 

+ / ( PI* SQRT (P I * S2 ) ) 

* * * IF SCATTERING FOIL IS PRESENT TWFOIL =0.5 AND SYMMETRIC TERM 

* * * IN IPERP AND IPAR IS DOUBLED 
RETURN 
END 

SUBROUTINE NORM (T, IPAR, IPERP, W,TAU, INORM, RMAX, QMAX, PMAX, AMAX, MAXY, 

+ MAXP) 

REAL INT ( 64 ) , IPAR (64) , IPERP (64) ,W(64) ,TAU(64) , PMAX 
REAL QMAX, AMAX, MAXY, MAXP, RMAX 
DO 100 NNW= 1 , 64 

IF (IPERP (NNW) .EQ. 0.0) THEN 
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noon nnnno 



TAU (NNW) =1 . 0 



ELSE 

TAU (NNW) = (IPERP (NNW) -I PAR (NNW) ) / (IPAR (NNW) +IPERP (NNW) ) 

TAU (NNW) =ABS (TAU (NNW) ) 

ENDIF 

C * * * * ADD IPAR & IPERP FOR UNPOLARIZED INTENSITY 

W (NNW) =IPAR (NNW) 4 IPERP (NNW) 

IF(IPAR(NNW) .GE.QMAX) THEN 
QMAX=IPAR (NNW) 

ENDIF 

IF (IPERP (NNW) .GE.PMAX) THEN 
PMAX=IPERP (NNW) 

ENDIF 

IF (W (NNW) . GE . RMAX) THEN 
RMAX=W (NNW) 

AMAX=T 

MAX Y= IPAR (NNW) 

MAXP= IPERP (NNW) 

ENDIF 

100 CONTINUE 

RETURN 
END 

SUBROUTINE BANDIN (KK, T, G, S, S2 , BMIN, BMAX, DPAR, DPERP ) 

EXTERNAL F 
DIMENSION P (5) 

DATA NN/3/ , FGI2/0 . 0/ 

***P(1) WILL CONTAIN VALUE OF DPERP IN F(B,P> CALLS, WHILE F ITSELF 
***IS VALUE OF DPAR 

* * *NN=NO . OF PTS . PER SUBINTERVAL FOR GAUSSIAN INTEGRATION 

* * * KK=NO . OF SUBINTERVALS IN BMIN TO BMAX, INTEG. RANGE 

* * * F IS NAME OF FUNCTION CALLED BY FGI : F=F(X,P) 

P (2) =T 
P (3)=G 
P (4 ) =S 
P (5) =S2 

DPAR=FGI (BMIN, BMAX, NN, KK, F, P, FGI2) 

DPERP=FGI 2 
RETURN 
END 

FUNCTION FGI (A, B, NN, K, F, P, FGI2) 

5.01.05 FGI - FORTRAN GAUSSIAN INTEGRATION ROUTINE 

5.01.05 DECK DATES 4/7/75 

5.01.05 TEST CASE 36 

DIMENSION V (165) , VI (90) , W (165) , W1 (90) , SUM (32) , P (5) , SUM2 (32) 
DIMENSION V2 (75) , W2 (75) 

EQUIVALENCE (V(l) , V2 (1) ) , (V<76) , VI (1) ) , 

4 (W(l) ,W2 (1) ) , (W (76) ,W1 (1) ) 

DATA V2 / 

4 -0.774596669, 0. , 0.774596669,-0.861136312,-0.339981044, 

4 .339981044, . 8 6 1 1 3631 2 , - . 9061 7 984 6, - . 5384 693 1 0, .0, 

+ 0.538469310, 0.906179846, -0.93246 95 14, -0.661209386, -0.2 38 619186, 

4 0.238619186, 0.661209386, 0.932469514,-0.949107912,-0.741531186, 

4 -0.405845151, 0. , 0.405845151, 0.741531186, 0.949107912, 
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» -0 . 960289856, 
+ 0.525532410, 

+ -0 . 613371 4 33 , 
+ 0.836031107, 
i- -0.4 33395394 , 
+ 0.865063367, 

+ -0.519096129, 
+ 0.730152001, 

+ -0.769902674, 
+ 0.367831499, 

DATA VI/ 

+ -0.984183055, 
+ -0.230458316, 
0.801578091, 

+ -0.827201315, 
+ 0.108054949, 

+ 0.928434884, 

+ -0.724417731, 
+ 0.201194094, 

+ 0.937273392, 

+ -0.755404408, 
+ 0.095012500, 
f 0.865631202, 

4 -0.964762256, 
^ -0.732182119, 
4 -0.331868602, 
4 0.144471962, 

4 0.587715757, 

4 0.896321156, 
DATA W2 / 

4 0.555555556, 

4 0.652145155, 

4 0.478628670, 

4 - 0.467913935, 

► 0.381830051, 

4 0.101228536, 

^ 0.313706646, 

4 0.260610696, 

4 0.180648161, 

4 0.269266719, 

4 0.149451349, 

4 0.233193765, 

4 0.186290211, 

4 0.160078329, 

4 0.233492537, 



-0 .796666477, 
0.796666477, 
-0 . 324253423, 
0.968160240,- 
-0 . 148874339, 
0 . 973906529, - 
-0.269543156, 
0 .887062600, 
-0.587317954, 
0.587317954, 



-0.525532410, 
0 . 960289856, - 

0 . 

0.973906529, - 
0.148874339, 
0 . 978228658, - 
0 . 

0 . 978228658, 
-0.367831499, 
0.769902674, 



-0 .183434642, 
0 . 968160240,- 
0.324253423, 
0.865063367,- 
0.433395394, 
0.887062600,- 
0.269543156, 
0 . 981560634, 
0.125233409, 
0.904117256, 



0.183434642, 
0. 836031107, 
0.613371433, 
0 . 679409568, 
0.679409568, 
0.730152006, 
0.519096129, 
-0.904117256, 
0.125233409, 
0.981560634/ 



-0 . 917598300, 
0 . 

0.917598400, 

-0.687292905, 

0.319112369, 

0.986283809,- 

-0.570972173, 

0. 394151347, 

0.987992518,- 

-0.617876244, 

0.281603551, 

0.944575023, 

-0 . 934906076, 

-0.663044267, 

-0.239287362, 

0.239287362, 

0. 663044267, 

0. 934906076, 

0 . 888888889, 

0 . 347854845, 
0.236926885, 
0.360761573, 
0.417959184, 
0.222381034, 
0.222381034, 

0 . 312347077, 
0.081274388, 
0.295524225, 

0 . 066671344, 
0.262804545, 
0. 125580369, 
0.203] 67427, 
0.203167427, 



-0.801578091, 
0.230458316, 
0 . 984183055, - 
-0.515248636, 
0.515248636, 

0 . 987992518, - 
-0.394151347, 
0.570972173, 

0 . 989400935, - 

-0.458016778, 

0 .458016778, 

0.989400935,- 

-0.896321156, 

-0.587715757, 

-0.144471962, 

0.331868602, 

0.732182119, 

0.964762256, 

0.555555556, 

0.236926885, 

0.171324492, 

0 .171324492, 
0.381830051, 
0.313706646, 
0.101228536, 

0 . 330239355, 

0 . 066671344, 
0.295524225, 

0 . 055668567, 

0 . 272925087, 

0 . 055668567, 
0.233492537, 

0 .160078329, 



-0.642349339, 
0.448492751, 
0.986283809,- 
-0 . 319112369, 
0.687292905, 
0.937273392,- 
-0.201194094, 
0.724417731, 

0 . 944575023, - 
-0.281603551, 
0 . 617876244, 

0 . 997263862,- 

-0.849367614, 

-0.506899909, 

-0.048307666, 

0 . 421351276, 

0.794483796, 

0 . 985611512, 

0.347854845, 

0.478628670, 

0.360761573, 

0.129484966, 

0.279705391, 

0.362683783, 

0.081274388, 

0 . 312347077 , 

0 .149451349, 

0 .269266719, 

0 . 125580369, 
0.262804545, 
0.047175336, 

0 .249147046, 

0 . 106939326, 



-0.448492751, 

0.642349339, 

0.928434860, 

-0.108054949, 

0.827201315, 

0.848206583, 

0 . 

0 . 848206583, 
0.865631202, 
-0 . 095012510, 
0.755404408, 
0.985611512, 
-0.794483796, 
-0.421351276, 
0.048307666, 
0.506899909, 
0.849367614, 
0.997263862/ 

0 . 652145155, 
0.568888889, 
0.467913935, 
0.279705391, 
0.129484966, 
0.362683783, 
0.180648161, 
0.260610696, 
0.219086363, 
0.219086363, 
0.186290211, 
0.233193765, 
0.106939326, 

0 . 2491 47046, 
0.047175336/ 



DATA Wl/ 

+ 0.040484005, 
4 0.226283180, 
+ 0.138873510, 
+ 0.121518571, 
+ 0.215263853, 
+ 0.080158087, 
+ 0.139570678, 
+ 0.198431485, 



0.092121500, 
0.232551553, 
0.092121500, 
0. 157203167, 
0.205198464, 
0.035119460, 
0. 166269206, 
0 . 186161000, 



0.138873510, 
0.226283180, 
0.040484005, 
0.185538397, 
0 .185538397, 
0.030753242, 
0 .186161000, 
0.166269206, 



0.178145981, 
0.207816048, 
0.035119460, 
0.205198464, 
0 . 157203167, 
0.070366047, 
0 .198431485, 
0 .139570678, 



0.207816048, 

0.178145981, 

0.080158087, 

0.215263853, 

0.121518571, 

0.107159220, 

0.202578242, 

0.107159220, 
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+ 0.070366047, 0.030753242, 0.027152459, 0.062253524, 0.095158512, 

+ 0.124628971, 0.149595989, 0.169156519, 0.182603415, 0.189450610, 

+ 0.189450610, 0.182603415, 0.169156519, 0.149595989, 0.124628971, 

+ 0.095158512, 0.062253524, 0.027152459, 0.007018610, 0.016274395, 

+ 0.025392065, 0.034273863, 0.042835898, 0.050998059, 0.058684093, 

+ 0.065822223, 0.072345794, 0.078193896, 0.083311924, 0.087652093, 

+ 0.091173879, 0.093844399, 0.095638720, 0.096540089, 0.096540089, 

+ 0.095638720, 0.093844399, 0.091173879, 0.087652093, 0.083311924, 
+ 0.078193896, 0.072345794, 0.065822223, 0.058684093, 0.050998059, 
+ 0.042835898, 0.034273863, 0.025392065, 0.016274395, 0.007018610/ 

N=MIN0 (MAXO (NN, 3) , 17) 

M=N 

IF (N.GE . 17) M=32 
NADD = 0 

IF (N.LE.3) GO TO 40 
IN = N-l 
DO 30 1=3, IN 
NADD = NADD + I 
30 CONTINUE 

40 DO 5 1=1, M 

SUM (I) =0 . 0 
SUM2 ( I ) =0 . 0 
5 CONTINUE 

11= (B-A) /FLOAT (K) 

H2=H/2 . 

AA=A+U2 

DO 20 L= 1 , K 

DO 10 1=1, M 

11 = NADD + I 

X = M2* V (II) + AA 

SUM ( 1 ) =SUM ( 1 ) +F (X, P) 

SUM2 (I) =SUM2 (I) +P (1) 

10 CONTINUE 

AA=AA+H 

20 CONTINUE 

SUMT=0 . 0 
SUMT2=0 . 0 
DO 25 1 = 1, M 
II = NADD + I 

SUMT = SUMT + W( II)* SUM (I) 

SUMT2 = SUMT2+W (II) * SUM2 ( I ) 

25 CONTINUE 

FGI=U2 * SUMT 
FGI2=U2* SUMT2 
RETURN 
END 

FUNCTION F (B, P) 

DIMENSION P (5 ) 

COMPLEX Q, RTQ, Z2, Z4 , AI 
Al= (0.0, 1.0) 

P (1) =0. 

T=P (2) 

G=P (3) 

S=P (4) 
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c 

c 

c 

c 

c 

c 

c 



c 



S2=P (5) 

Q=1 . 0/S2 -AI * 2 . 0 * B 
RTQ=CSQRT (Q) 

Z2=AI * RTQ/G 
Z4=Q/G 

CALL DIST2 (T, G, S, Q, B, S2 , RTQ, Z2 , Z4 , DPAR, DPERP) 

F=DPAR 
P (1) =DPERP 
RETURN 

end 

SUBROUTINE CERF (MO, A, W) 

****************************************************************** 

FORTRAN SUBROUTINE FOR COMPLEX ERROR FUNCTION 
****************************************************************** 

MO = MODE OF OPERATION 
A = ARGUMENT (COMPLEX NUMBER) 

W = FUNCTION (COMPLEX NUMBER) 



COMPLEX A, W 

DIMENSION AZ (2) , QZ (2) , SZ<2), EF(2) 

DIMENSION TS (2 ) , SM(2), TM<2), QF<2),. CD(18), CE(18) 

DATA CD (1) /0. 00000000000000 E00/, CD ( 2 ) /2 . 0 860585 601 347 6E-2/ 
DATA CD (3) /8. 298069404 95 687 E- 2/, CD ( 4 ) / 1 . 854 21 65332 6079E-1 / 
DATA CD(5) /3.27963479382361E-1/,CD(6) /5 . 12675279912828E-1 / 
DATA CD (7) /I . 45412958045105E-1/, CD (8) /I .0369506741 8297E00/ 
DATA CD (9) /I . 40378061255437E00/,CD(10) /I . 86891662214001E00/ 
DATA CD (11) /2. 4631483052392 9E 00/, CD (12) /3 . 2271 93837 37 352E00/ 
DATA CD (13) / 4 .21534348280013E00/,CD (14) /5 . 50 17 8 87 3 151 54 9E00/ 
DATA CD (15) /7 . 19258966683102E00/, CD (16) /9 . 45170208076408E00/ 
DATA CD (17) /I .25710718314784E+1/, CD (18) /I . 724 835372 16334E+1/ 
DATA CE (1) /8 . 1 572308 3324096E-2/ , CE ( 2 ) /I . 592852852534 37E-1 / 
DATA CE (3) /I . 4858162561449 9E-1/, CE ( 4 ) / 1 . 332 1 96708362 4 5E- 1 / 
DATA CE(5) /I . 15690392878957 E-1/,CE(6) /9 . 78580959447535E-2/ 
DATA CE (7) /8 .059088342 97624E-2/,CE (8) /6 . 4 0204 538 609872E-2/ 
DATA CE (9) / 4 .8144 524 27 67885E-2/ , CE ( 1 0 ) /3 . 3354 065847 32 95E-2/ 
DATA CE (11) /2. 0554 8099470 193E-2/,CE (12) /I . 07847403887506E-2/ 
DATA CE (13) /4 .55634 8922142l9E-3/,CE (14) /I .4 39844581 38925E-3/ 
DATA CE (15) /3 . 07 0561 39834 1 7 IE-4 / , CE (16) /3 . 78 1 5654 1 1 6854 IE-5/ 
DATA CE (17) /2 . 05 17 350961 61 2 IE- 6/, CE (18) / 2 . 63564 82368274 7E-8/ 



AZ (1) =REAL(A) 

AZ (2) =AIMAG (A) 

001 ZS=AZ (1) *AZ (1) +AZ (2) *AZ (2) 

SZ ( 1 ) =AZ ( 1 ) *AZ ( 1 ) -AZ (2) *AZ (2) 
SZ (2) =2 . 0 * AZ (1) *AZ ( 2 ) 

QZ ( 1 ) =+AZ ( 1 ) 

QZ (2) =+AZ (2) 

SN=+1 . 0 

IF (AZ (1) ) 002, 003, 003 

002 QZ ( 1 ) =-AZ ( 1 ) 

QZ (2) =-AZ (2) 

SN=-1 .0 

003 IF (ZS-1 .0) 014, 014, 004 

004 IF (ZS-38 . 0) 005, 006, 006 
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ODD 



005 IF (SZ (1) +0. 0 64 *SZ (2) *SZ (2) ) 014, 014, 012 
C 

C ASYMPTOTIC EXPANSION 

C 

006 SM(1)=0.0 
SM (2) =0 . 0 

QF (1) =4SZ (1) / (ZS*ZS) 

QF (2) =-SZ (2) / (ZS*ZS) 

QM=0 .5641 8 958354 77 56* EXP (-SZ (1) ) 

TS (1) =QM*COS (-SZ (2) ) 

TS (2) =QM*SIN (-SZ (2) ) 

TM ( 1 ) =+ (TS ( 1 ) *QZ ( 1 ) +TS (2 ) *QZ<2) ) /ZS 
TM(2) =- (TS (1) *QZ (2) -TS (2) *QZ (1 ) ) /ZS 
PM=-0 . 5 
GO TO 008 

007 PM=PM+ 1 . 0 

TS (1) =TM(1) *QF (1) -TM<2) *QF (2) 

TS (2) =TM(1) *QF (2) 4TM (2) *QF (1) 

TM (1 ) =-PM*TS ( 1 ) 

TM (2) =-PM *TS (2) 

IF (ABS (SM ( 1 ) ) + ABS (TM ( 1 ) ) . NE . ABS ( SM ( 1 ) ) ) GO TO 008 
IF (ABS (SM (2) ) +ABS (TM (2) ) . EQ . ABS (SM(2) ) ) GO TO 009 

008 SM ( 1 ) =SM ( 1 ) + TM ( 1 ) 

SM (2) =SM (2 ) t TM (2) 

IF (PM . LT . 36 . 5 ) GO TO 007 

009 IF (QZ(1) .LT.0.01) GO TO 019 
EF ( 1 ) =SM ( 1 ) 

EF (2) =SM (2) 

010 IF (MO . EQ . 0) GO TO Oil 
W=CMP LX (EF (1) , EF (2) ) 

IF ( S N . EQ .1.0) RETURN 
W=CMPLX (2 . 0-EF (1) , -EF (2) ) 

RETURN 

011 EF (1) =SN* (1 . 0-EF (1 ) ) 

EF (2) =-SN*F,F (2) 

W=CMP LX (EF (1) , EF (2) ) 

RETURN 

RATIONAL FUNCTION APPROXIMATION 

012 SM ( 1 ) =0 . 0 
SM (2 ) =0 . 0 

QM=0 .564 18958354 7756* EXP (-SZ (1) ) 

TS ( 1 ) =QM*COS (-SZ (2 ) ) 

TS (2) =QM*SIN (-SZ (2) ) 

QF ( 1 ) =TS ( 1 ) *QZ(1)-TS(2) *QZ (2) 

QF (2)=TS(1) *QZ(2) +TS (2) *QZ(1) 

DO 013 1=1,18 
TS (1) =SZ (1) 4CD (I) 

TS (2) =SZ (2) 

SS=TS (1) *TS (1) 4TS (2) *TS (2) 

TM ( 1 ) =+CE ( I ) * TS (1) /SS 
TM ( 2 ) =-CE ( I) *TS (2) /SS 
SM ( 1 ) =SM ( 1 ) I TM ( 1 ) 
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SM ( 2 ) =SM ( 2 ) +TM ( 2 ) 

013 CONTINUE 

EF ( 1 ) =QF (1) * SM ( 1 ) -QF (2) *SM(2) 

EF ( 2 ) =QF (1) *SM(2) +QF (2) *SM(1) 

GO TO 010 

TAYLOR SERIES EXPANSION 

014 TM(1)=1 . 12837 91 67095513 *QZ (1) 

TM(2)=1. 1283791 67095513*QZ (2) 

SM ( 1) =TM (1) 

SM ( 2 ) =TM ( 2 ) 

PM=0 . 0 

015 PM=PM+ 1 . 0 
DM=2 . 0*PM+1 .0 

TS ( 1 ) =TM ( 1 ) *SZ (1) -TM(2) *SZ (2) 

TS ( 2 ) =TM ( 1 ) * SZ ( 2 ) +TM ( 2 ) * SZ ( 1 ) 

TM< 1)=-TS (1) /PM 
TM ( 2 ) =-TS (2) /PM 
TS ( 1 ) =TM (1) /DM 
TS (2) =TM (2) /DM 

IF (ABS (SM (1 ) ) 4 ABS (TS ( 1) ) -NE . ABS (SM (1) ) ) GO TO 016 
IF (ABS (SM (2) ) +ABS (TS (2) ) . EQ . ABS (SM (2) ) ) GO TO 017 

016 SM ( 1 ) =SM ( 1 ) +TS ( 1 ) 

SM (2) =SM (2) +TS (2) 

GO TO 015 

017 I F (MO . NE . 0 ) GO TO 018 
W=CMPLX ( SN * SM ( 1 ) , SN * SM (2 ) ) 

RETURN 

018 EF ( 1) =1 . 0-SN*SM (1) 

EF (2) =-SN*SM (2) 

W=CMPLX (EF (1) , EF (2) ) 

RETURN 

MODIFIED ASYMPTOTIC EXPANSION 

019 SN=-SN 
GO TO 017 
END 

"Z 
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